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ABSTRACT 

Many of the carbon-enhanced metal-poor (CEMP) stars that we observe in the Galactic halo are found in binary systems and show 
enhanced abundances of elements produced by the slow neutron-capture process (j-elements). The origin of the peculiar chemical 
abundances of these CEMP-i stars is believed to be accretion in the past of enriched material from a primary star in the asymptotic 
giant branch (AGB) phase of its evolution. We investigate the mechanism of mass transfer and the process of nucleosynthesis in 
low-metallicity AGB stars by modelling the binary systems in which the observed CEMP-i- stars were formed. For this purpose we 
compare a sample of 67 CEMP-i stars with a grid of binary stars generated by our binary evolution and nucleosynthesis model. We 
classify our sample CEMP-i stars in three groups based on the observed abundance of europium. In CEMP -^/r stars the europium-to- 
iron ratio is more than ten times higher than in the Sun, whereas it is lower than this threshold in CEMP-.s/nr stars. No measurement 
of europium is currently available for CEMP-j/ur stars. On average our models reproduce well the abundances observed in CEMP- 
s/nr stars, whereas in CEMP-s/r stars and CEMP-j/ur stars the abundances of the light-s elements (strontium, yttrium, zirconium) 
are systematically overpredicted by our models and in CEMP-s/r stars the abundances of the heavy-s elements (barium, lanthanum) 
are underestimated. In all stars our modelled abundances of sodium overestimate the observations. This discrepancy is reduced only in 
models that underestimate the abundances of most of the ^-elements. Furthermore, the abundance of lead is underpredicted in most of 
our model stars, independent of the metallicity. These results point to the limitations of our AGB nucleosynthesis model, particularly 
in the predictions of the element-to-element ratios. In our models CEMP-j stars are typically formed in wide systems with periods 
above 10,000 days, while most of the observed CEMP-i stars are found in relatively close orbits with periods below 5,000 days. This 
evidence suggests that either the sample of CEMP-i binary stars with known orbital parameters is biased towards short periods, or 
that our wind mass-transfer model requires more efficient accretion in close orbits. 
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1. Introduction 

The very metal-poor stars observed in the Galactic halo are 
of low mass and exhibit abundances of iron of approximately 
[Fe/H] < -2.0. Very metal-poor stars carry the fingerprints of 
the early stages of evolution of the Milky Way and therefore 
have been extensively studied by different surveys in the past 
two decades, with particular focus on chemically peculiar stars 
(e.g.:|Beers et al.|1992[|Christlieb et al.|2001[|Frebel et al.|2006| 

|Yanny et al.||2009| l. Among these, carbon-enhanced metal-poor 

(CEMP) stars are very metal-poor stars enriched in carbon. The 
observed fraction of CEMP stars in the halo increases with in¬ 
creasing galactic latitude and with decreasing metallicity, and 
varies between approximately 9% at [Fe/H] < -2 and about 


25% at [Fe/H] < -3 (jCohen et al. 

2005 

Marsteller et al.|20051 

Frebel et al.|2006| Lucatello et al. 

2006 

Lee et al.||2013 

Yong 

et al.|2013|l. 


In the literature CEMP stars are generally defined by the 
carbon excess [C/Fe] > 1.0, and are classified in groups ac¬ 
cording to the observed abundances of barium and europium, 
two heavy elements produced by the slow (s-) and the rapid (r-) 


neutron-capture process, respectively. The exact definitions vary 
between different authors (e.g. Beers & Christlieb|2005 |Jonsell 


|et al.|2006l|Aoki et al.|20()7l|Masseron et al.|2010] l, and in this 

work we adopt the following classification scheme. 


/. CEMP-s stars are CEMP stars that satisfy the criteria 
[Ba/Fe] > 0.5 and [Ba/Eu] > 0. 
a. CEMP-s/r stars are CEMP-s stars enriched in europium, i.e. 
[Eu/Fe] > 1. 

Hi. CEMP-r stars have [Eu/Fe] > 1 and [Ba/Eu] < 0. 
iv. CEMP-no stars do not exhibit enhanced abundance of bar¬ 
ium, i.e. [Ba/Fe] < 0.5. 


It has been suggested that CEMP stars owe their abundances 
to mass transfer of carbon-rich material in the past from a 
thermally-pulsing asymptotic giant branch (TP-AGB) primary 
star that today is an unseen white dwarf (e.g. Wallerstein & 


Knap p|1998[ Preston & Sneden|2001 Beers & Christlieb|2005 
Ryan et al.|2005| ). 

In this work we focus on CEMP-s stars, for which there are 
stronger arguments in favour of the binary mass-transfer sce¬ 
nario. Qualitatively, in this scenario the primary star evolves to 
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the AGB phase, produces carbon and heavy elements and sub¬ 
sequently transfers part of this material by wind mass transfer to 
the low-mass companion star, the star observed today. A quanti¬ 
tative model of this process depends on many aspects of stellar 
evolution, AGB nucleosynthesis and binary interaction that are 
not well understood. In AGB stars many uncertainties are related 
to the physics of mixing, which determines the stellar structure 
and chemical composition and which in turn influence the radius 
of the star, and hence its luminosity and mass loss rate (|Herwig 


|2005 1 l^nstantino et al. |20141 IFishlock et al. |20 1 4| l . The material 
that is expelled by the AGB star carries away angular momen¬ 
tum from the system and its fate depends on the mass transfer 
mechanism. Therefore the study of CEMP-i stars provides us 
with a powerful tool to improve our understanding of AGB nu¬ 
cleosynthesis at low metallicity and of the mass-transfer process 
in binary systems. 

In our previous paper ( Abate et al.|20I5a hereinafter Paper I) 
we analyse a sample of 15 CEMP binary stars with known orbital 
periods: through the comparison with the observed abundances 
while matching the measured periods we put new constraints on 
our models of binary stellar evolution and AGB nucleosynthe¬ 
sis. In most of the systems a combination of large mass accre¬ 
tion and efficient angular momentum loss is necessary to repro¬ 
duce at the same time the observed chemical abundances and 
orbital periods. About half of the stars of the sample are not ac¬ 
curately reproduced by any of our models, regardless of the as¬ 
sumptions made about the dynamical evolution of the systems, 
and this points to the limitations in our AGB nucleosynthesis 
model, particularly on the maximum abundances of the heavy 
i-process elements and on the element-to-element ratios that are 
produced. 

In this work we extend the analysis of Paper I to a larger sam¬ 
ple of 67 CEMP-s stars that includes systems without informa¬ 
tion about the orbital period. We model a grid of about 400,000 
binary stars with different masses and separations and we com¬ 
pare the modelled abundances with the observations. Eor each 
star in the sample we determine the model star that best matches 
the observed abundances, with the same procedure of mini¬ 
mization followed in Paper I. In every star for each observed el¬ 
ement we compute the residual as the difference between the ob¬ 
served and the modelled abundance and subsequently we anal¬ 
yse the distributions of the residuals for every element individu¬ 
ally. The purpose of this analysis is to investigate if statistically 
our model reproduces the observations, even though discrepan¬ 
cies may occur for some elements in individual stars, and thus to 
test our models with a study complementary to that of Paper I. 

The paper is organised as follows. In Sect.j^we describe our 
observational sample, and in Sect. we present our model of 
binary stellar evolution and nucleosynthesis. In Sect. we dis¬ 
cuss the set of initial abundances adopted in our model. In Sect. 
I^we present the results of the analysis of our sample CEMP-s 
stars, we compare the modelled and observed abundances of ev¬ 
ery element individually and we discuss the distribution of the 
residuals. In Sect. we investigate the confidence limits on the 
initial parameters of our best fits and we study the distribution of 
these parameters. In Sect. we discuss our results while Sect.j^ 
concludes. 


2. Data sample 



CEMP-s/nr 
CEMP-s/ur 
CEMP-s/r 
CEMP-no 
no Ba, Eu 
C-normal 


X 

□ 


Fig.l. [C/Fe] vs [Fe/H] in the stars of our observed sample. 
Grey plus signs (-t) are carbon-normal, very metal-poor stars; 
CEMP-i stars are classified in three groups according to the 
abundance of europium, as described in the text: CEMP-s/nr 
stars (•), CEMP-i/wr stars (x) and CEMP-i/r stars (□). Purple 
triangles (a) and grey dots (.) are CEMP-no stars and CEMP 
stars without measured abundances of heavy elements, respec¬ 
tively. The dotted line indicates [C/Fe] = 1. 


with observed abundances of carbon and barium and we ignore 
stars with only upper or lower limits. In some stars measure¬ 
ments of element abundances or stellar parameters are available 
from multiple sources. In most cases the measurements are con¬ 
sistent within the observational uncertainties and for the purpose 
of our study we use the arithmetic mean of the logarithm of the 
observed abundances and we adopt the largest observed error as 
the uncertainty on the measure. In case two measurements differ 
by more than the declared observational uncertainty and there is 
no obvious criterium to prefer one value, we compute the aver¬ 
age and we adopt as the uncertainty half the difference between 
the two values. To this sample we add four CEMP stars studied 
by Masseron et al. ( 2010| l with metallicity in the above range 
that were not present in the SAGA database. 

This selection leaves us with a sample of 378 very metal- 
poor stars, 67 of which are classified as CEMP-i stars, 8 as 
CEMP-no stars, 46 are CEMP stars with no information about 
the abundances of neutron-capture elements, and the remaining 
257 are carbon-normal very metal-poor stars. In Fig. we show 
the carbon abundances of the stars in our observed sample as a 
function of [Fe/H]. CEMP-no stars are indicated as grey trian¬ 
gles. The dotted line represents the threshold carbon abundance 
[C/Fe] = 1 above which the stars are defined CEMP stars. We 
classify CEMP-s stars in three groups based on the abundance 
of europium: 


• In CEMP-s/r stars the abundance of europium is enhanced, 
[Eu/Fe] > 1 (open squares). 

• In CEMP-s/nr stars the europium abundance is [Eu/Fe] < 1 
(filled circles). 

• In CEMP-s/ Mr stars the abundance of europium has not been 
determined (crosses), for example because the spectra have 
low signal-to-noise ratios, or the europium lines are blended. 


Our database of observed very metal-poor stars is based on 580 


stars catalogued in the SAGA observational database (Suda et al. 
2008 201 1[ last updated in January 2015) with iron abundance 


-2.8 < [Fe/H] < -1.8. Among these objects we select the stars 


Table |A.1| summarises the surface gravities, temperatures, and 
abundances of iron, carbon, barium and europium (when avail¬ 
able) of the 67 CEMP-s stars in our observed sample. For three 
stars in which barium is not observed and lanthanum is en- 
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hanced (HD13826, HD198269 and HD201626) the abundance 
of lanthanum is listed. The minimum uncertainty assumed in the 
chemical abundances and surface gravities is 0.1 dex. We adopt 
an uncertainty of 100 K in effective temperature unless differ¬ 
ently stated. 


3. Models of binary evolution and nucleosynthesis 


As in Paper I, in this study we use the code binary_c/nucsyn 
that couples algorithms to compute the evolution of stars in bi¬ 
nary systems with a model of stellar nucleosynthesis. The details 
of our code and the prescriptions used for the binary stellar evo¬ 


lution and nucleosynthesis are extensively discussed by Izzard 
et al. ( 2004| 2006 2009| l. Our default input physics is the same 
as tn Paper I, and we refer to this for a more complete descrip¬ 
tion. In this section we summarise some important parameters 
adopted in our model, related to the wind mass-transfer process 
and the nucleosynthesis in the AGB phase (Sect. 3.1 and 3.2 1 , 


and we describe the basic characteristics of our grid of model 
stars (Sect.[X3|l. 


3.1. Wind-accretion rate and anguiar momentum ioss 

In this work we compare the results derived with three differ¬ 
ent model sets based on different assumptions about the wind- 
accretion rate and the mechanism of angular momentum loss. 
The model sets are listed in Table [T] Model sets A and B are 
the same as in Paper I. In our default model set A we calculate 
the wind-accretion rate according to the prescription for wind 
Roche-lobe overflow (WRLOF) with a dependence on the mass 
ratio, as presented by |Abate et al.| ( |2013| Eq. 9). We compute the 
angular momentum carried away by the wind material assuming 
a spherically symmetric wind (Eq. 4 of Abate et al.|2013 1. The 
results of Paper I show that a mechanism of wind accretion that 
is much more efficient at relatively short separations is neces¬ 
sary to reproduce the abundances observed in two thirds of the 
analysed CEMP-i stars. To take this into account, in model set 
B we adopt the Bondi-Hoyle-Lyttleton (BHL) prescription as in 
Eq. ( 6 ) of Boffin & Jorissen ( |1988 1 with asHL = 10 (instead of 
1 < aBHL 5 2 ) to simulate a very efficient mass-transfer process. 
Model set B also adopts a prescription of efficient angular mo¬ 
mentum loss, in which the material lost from the binary system 
carries away a multiple y = 2 of the average specific orbital an¬ 
gular momentum (as in Eq. 2 of Izzard et al., 2010, and Eq. 10 
of |Abate et al.|201^ . Model set C combines the WRLOF accre¬ 
tion rate as in model set A with the efficient angular momentum 
loss as in model set B. In the analysis performed by Abate et al. 
(|2013|l this is the model set that predicts the largest fraction of 
CEMP stars. 


3.2. AGB nucieosynthesis 


We refer to the review papers by Busso et al. ( 1999[ ) and |Herv^ 
( |2005| l for a description of the nucleosynthesis process in the 
interior of AGB stars, the role of the ^^C-pocket at the top of 
the intershell region as a source of free neutrons for the produc¬ 
tion of slow neutron-capture elements (s-elements) in low-mass 
AGB stars, and the effect of the third dredge-up process (TDU 
hereinafter) that mixes the products of i nternal nucleosy nthesis 


Karakas 


(120101 a i^c- 


to the stellar surface. In the models of 
pocket is created by including a partial mixing zone (PMZ) at 
the deepest extent of each TDU, where protons are mixed in the 
intershell region and subsequently captured by the nuclei to 


Table 1. Models of the wind-accretion efficiency and angular 
momentum loss used in this study. 


Model set 

wind accretion 

angular momentum 


efficiency 

loss 

A 

WRLOE 

spherically symmetric wind 

B 

BHL, ttBHL - 10 

AJ/J = 7 (AM/M) ,7 = 2 

C 

WRLOE 

AJ/J = 7 (AM/M) ,7 = 2 


Notes. Our WRLOF model is calculated with Eq. (9) of|Abate et al. 

1 2013 1 . The BHL model is computed with Eq. (6) of |Boffin & Jorissm 
1 1988 1 . Angular momentum loss due to wind ejection is computed alter¬ 
natively with Eq. (4) of |Abate et al.|(|2013) for a spherically symmetric 
wind or with y = 2 in Eq. (2) of [Izzard et al.| l |2010) and Eq. (10) of 
[Abate et al.| ( |2013l . 


form a layer rich in The mass of the PMZ is a free param¬ 
eter in the models, and |Lugaro et al. ( 2012[ ) study the effects of 
different masses of the PMZ on the surface abundances in AGB 
stars with different initial masses. 

We describe in Paper I how in our models the amount of 
material dredged-up from the intershell region is determined in 
order to reproduce the evolution predicted in the detailed mod¬ 
els of Karakas ( |2010 l and |Lugaro et al. ( 2012| l. The chemical 
composition of the intershell region is saved in a table as a func¬ 
tion of three parameters: the mass of the star at the beginning 
of the TP-AGB phase, the thermal-pulse number, and the mass 
of the PMZ, MpMz- The surface abundances of an AGB star of 
mass M, evolve in time and are recalculated at every TDU, when 
material with the chemical composition of the intershell region 
taken from our table is mixed into the convective envelope. 


3.3. Grid of modeis 


Our simulations are based on the same grid of models as in 
Paper I, with N binary-evolution models distributed in the Mi - 
M 2 - logio a - Mpuz parameter space, where Mi ,2 are the initial 
masses of the primary and secondary stars, respectively, a is the 
initial separation of the system, and Mpmz is the mass of the par¬ 
tial mixing zone of any star of the binary system that becomes an 
AGB star. The grid resolution N - Nmi xNmiXNoXNpmz^ where 
we choose Ami = 34, Am 2 = 28, Na = 30, and Apmz = 10. The 
initial parameters are chosen as follows: 


Ml varies in the range [0.9,6.0] Mq. The grid spacing is 
AMi i = 0.1 Mq up to 3 Mg, and AMi i = 0.25 Mg otherwise. 
M 2 is equally spaced in the range [0.2,0.9] Mg, and by defi¬ 
nition M 2 ,i < Mi l. 

a varies between 10^ Rg and 10^ Rg. The distribution of sep¬ 
arations is flat in logjg a. In the mass range considered here, 
at shorter separations the evolution of primary stars is in¬ 
terrupted before the AGB, because the systems undergo a 
common-envelope phase, whereas stars at wider separation 
do not interact in our models. The eccentricity is always zero. 
When Ml 1 < 3 Mg we adopt 10 different values for Mpmz. 
namely 0,’l0-\ 2 x lO^^, 5 x 10-\ 6.66 x lO-'*, 10-^ 1.5 x 
10“^, 2 X 10“^, 3 X 10“^, and 4 x 10“^ Mg. Mpmz is zero 
otherwise, in accordance with the detailed models of AGB 
nucleosynthesis of Karakas (2010|l. 


All stars in our grid are formed in binary systems. In reality 
single stars or binaries in wider orbits may exist, and carbon-rich 
very metal-poor AGB stars have recently been observed ( |Boyer 
jet al. 112015] ). To determine whether the abundances observed in 
the stars of our sample can be reproduced by a single-star model. 
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Fig. 2. Average abundances observed in our sample of carbon-normal stars ([C/Fe] < 1). Circles and plus signs represent elements 
observed in more, or fewer, than six stars, respectively. The solid and dashed bars enclose, respectively, 68.3% and 100% of the 
observations. The combined abundance of Ch-N is shown in the right panel. The solid line shows the initial-abundance mix adopted 
in our models. 


we add to our grid 200 synthetic single stars with masses uni¬ 
formly spaced in the range [0.4, 1.0] Mq, and Mpmz chosen as 
above. 

We assume that the accreted material is instantaneously 
mixed throughout the entire secondary star, as in the default 
model of Paper I. This approximation mimics the effect of ef- 
hcient non-convective mixing processes, such as thermohaline 
mixing, and is reasonable in low-mass stars ([Stancliffe et al.| 
[ 2 ^ . 


In our mod el we assume the sa me metallicity as in the de¬ 
tailed models of Lugaro et al. (20121, Z = 10 that corresponds 
approximately to [Fe/H] -2.3. As initial chemical composi¬ 
tion of the stars in our grid we adopt the abundances predicted in 
the one-zone Galactic chemical-evolution model of |Kobayashi| 
|et al.| ( |^ll| l at [Fe/H] = -2.3 that includes the contrihutions 
of core-collapse and type la supernovae and AGB stars. Their 
results extend up to ^®Ge and for heavier isotopes we adopt the 
solar abundance distribution derived by [Asplund et al. ( |2009| l 
scaled to metallicity Z = 10“^. Different assumptions on the ini¬ 
tial composition of the stars negligibly affect AGB nucleosyn¬ 
thesis, as discussed by Lugaro et al. ( |2012 1 . In the next section 
we discuss how representative our set of initial abundances is 
of the chemical composition of the very metal-poor stars in our 
sample. 


4. Comparison of the model initial abundances with 
carbon-normal metal-poor stars 

In this section we focus on the 257 “carbon-normal” very metal- 
poor stars in our sample, i.e. with [C/Fe] < 1, called C-normal 
stars hereinafter. C-normal stars do not exhibit evidence of du¬ 
plicity and are not expected to have changed their initial surface 
composition. Consequently, their abundances are expected to 
follow the distribution predicted by Galactic chemical-evolution 
models at metallicity Z a; 10 Some abundance variations 
may be introduced by internal mixing processes, such as rota¬ 
tional mixing, gravitational settling of heavy nuclei or thermo¬ 
haline mixing. The only signihcant change in the abundances 
of stars with small surface gravity, log[o(g/cms“^) < 4, is be¬ 
cause of the first dredge-up that reduces the carbon-to-nitrogen 
ratio but leaves essentially unaltered the abundances of the other 
elements. Our purpose is to compare the set of initial abun¬ 
dances adopted in our models with the abundances observed in 


C-normal stars and find discrepancies that are relevant to our 
study of the chemical composition of CEMP-i stars. 

In Fig.|^we show the initial abundances adopted in our mod¬ 
els (solid line) and the abundances observed in the C-normal 
stars. For every observed element we compute the mean of the 
logarithm abundance relative to iron, [El/Fe] (filled circles and 
plus signs). Only elements observed in more than one star are 
shown. Elements observed in fewer than seven stars are plotted 
with dashed bars, which connect the minimum and the maximum 
observed value. Eor each element observed in at least seven stars 
we calculate the median of the values of [El/Ee] and we select 
the group of approximately 68% stars, half of which have [El/Ee] 
larger than the median and the other half have smaller values. 
The solid bars in Fig. [^indicate the interval of [El/Fe] that en¬ 
closes the stars in this group. In the absence of signihcant nu- 
cleosythesis in C-normal stars, we expect that: (1) the mean and 
the median of the observed abundances coincide and are equal 
to the initial value predicted by the models of Galactic chemical 
evolution; (2) the deviations from the mean have a Gaussian dis¬ 
tribution; (3) the solid bars are symmetrical with respect to the 
mean and correspond to the standard deviation in the observed 
abundances. However, Fig. shows that these three statements 
are not always true, and several discrepancies occur between the 
predictions of the Galactic chemical evolution model and the ob¬ 
servations. 

We briehy analyse the discrepancies that are most relevant 
for our purposes. The differences between the model abundances 
and the observations of carbon and i-elements (e.g. Sr, Y, Zr, Ba, 
La, Ce, Pb) are on average within a factor of two (0.3 dex), al¬ 
though in some cases the observations have a large spread, for 
example carbon and barium (0.8 dex), cerium (0.9 dex) and lead 
(1.1 dex). These uncertainties are expected to have a small ef¬ 
fect on our study of CEMP-i stars because typically in low-mass 
AGB stars (M, < 3Mq) the amount of carbon increases by more 
than two orders of magnitude and the abundances of s-elements 
by more than a factor of ten. The average abundance of nitro¬ 
gen relative to iron is 1.5 dex larger than the initial abundance 
assumed in our model, with a spread of 1 dex. This discrep¬ 
ancy may alter the results of our best-fitting models. In stellar 
nucleosynthesis, nitrogen is made in the CN cycle when pro¬ 
tons are mixed into regions of the star where carbon is abundant. 
For example, the mixing of protons can occur at the bottom of 
the convective envelope of red giants and AGB stars during hrst 
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dredge-up and TDU, respectively. The increase in nitrogen abun¬ 
dance depends on the extent of the mixing, which is very uncer¬ 
tain (e.g.: Charbonnel et al.|1998[ Boothroyd & Sackmann|199^ 


[Karakas et al.|2010[ ). However, the total amount of carbon plus 

nitrogen is conserved in the CN cycle and therefore when both 
the elements are observed it is convenient to compare the pre¬ 
dictions of our models with the combined abundance C-hN. In 
C-normal stars the average Ch-N (right panel of Fig. ^ is un¬ 
derestimated by 0.5 dex (approximately a factor of three). This 
offset should have a small effect on our fits of CEMP-s stars, be¬ 
cause the abundance of Ch-N increases by more than a factor of 
a hundred in AGB stars. 


Oxygen is underpredicted on average by a factor of three. 
The amount of oxygen increases by less than one order of mag¬ 
nitude during the AGB phase, hence an initial offset can sub¬ 
stantially modify our results in the study for CEMP-s stars. The 
observed abundances of sodium are underpredicted on average 
by 0.3 dex and are spread over 0.7 dex. A better model of the 
initial abundance of sodium is desirable because the amount of 
sodium produced in AGB nucleosynthesis depends strongly on 
the stellar mass and on the mass of the PMZ, and a reliable fit of 
its abundance can, in principle, provide constraints on these two 
parameters. 


The elements between atomic numbers 13 (aluminium) and 
30 (zinc) are not always consistent with our initial assumptions. 
These elements are not modified by AGB nucleosynthesis, with 
the exception of aluminium that is partly produced in massive 
AGB stars, hence the offset observed in C-normal stars is ex¬ 
pected to be an issue also in CEMP-s stars. 


The elements with atomic numbers 62-71 and 75-80 are un¬ 
derestimated by 0.4-1 dex, in many cases with large spreads in 
the observed abundances. In the solar system, more than 60% 
of the total abundance of these elements is produced by the r- 
process, for example 69% of samarium, 94% of europium and 
86 % of gadolinium ( |Arlandini et al.|1999[ Bisterzo et al.|20lT i. 
Eor brevity, these elements are hereinafter referred to as “r- 
elements”, although some fraction of their total amount is pro¬ 
duced by the s-process. However, during AGB nucleosynthesis 
their abundances generally increase by less than 1-1.5 dex, and 
therefore the initial offset may affect our analysis of the CEMP-s 
stars. 


5. Analysis of the abundances in carbon-enhanced 
metal-poor stars 

We compare our binary evolution and nucleosynthesis models 
with the abundances of the 67 CEMP-s stars in our observational 
sample. We also include 14 of the 15 binary stars with known 
orbital periods studied in Paper I (those with -2.8 < [Ee/H] < 
1 .8) but in the present work, for comparison, we focus only on 
the chemical abundances, and we ignore the constraints on the 
period of the systems. In general we find initial parameters of 
the best-fit model different from Paper I. These differences give 
an estimate of the uncertainties in our present results caused by 
the lack of information on the orbital periods. 

Eor each star in our sample we determine the model in our 
grid that best reproduces the observed chemical abundances with 
the same procedure as in Paper I. Initially, to constrain the evolu¬ 
tionary stage, we select model stars that reproduce the measured 
surface gravity within the observational uncertainty, cTg, at an 
age 10 < f < 13.7 Gyr. Subsequently, for the stars that pass this 


selection, we compute the of each model from 




Z 


(^/,obs ^/,mod) 


cr 


(1) 


/, obs 


In Eq. Q every element i has observed abundance A, obs = 12 H- 
log[g(A,;obs/NH), where A;,obs and Ah are the number densities 
of i and hydrogen, respectively, cr, obs is the observational error 
associated with A,_obs and A, mod is the abundance predicted by 
the model. The minimum value of . , determines the best 
model. In this procedure we do not include any constraint on the 
effective temperature, Tetr, because Toff depends strongly on the 
observed metallicity which varies by up to a factor of three in the 
stars of our sample compared to the value adopted in our model, 
Z = 10“^. All other parameters being equal, if the observed 
metallicity is lower (higher) than in our model we expect to find 
a model Teff lower (higher) than the observed. 

As we noted in Paper I, because in our study we adopt a 
fixed metallicity for all our systems we ignore the scatter in ob¬ 
served [Ee/H] values. Qualitatively, a decrease in metallicity has 
two effects. Eirst, the abundance relative to iron of the elements 
produced in AGB nucleosynthesis increases, because the ini¬ 
tial amount of iron is lower, and iron is not produced in AGB 
stars. Second, because the production of neutrons is primary in 
AGB stars, the smaller the abundance of iron, the larger is the 
neutron-to-iron ratio, and consequently the abundances of the 
more neutron-rich elements (e.g. barium, lead) are enhanced. 
Opposite effects are caused by increasing metallicity. We re¬ 
stricted our sample to a 0.5 dex range around [Ee/H] = -2.3 to 
have a sufficient number of stars for our analysis. It is important 
to consider that our approximation probably introduces bigger 
errors the larger the difference in [Ee/H] between observations 
and our model. An attempt to quantify this error is discussed in 
Sect.im 

We focus on the elements produced by nucleosynthesis in 
AGB stars. Therefore in Eq. Q we take into account C, N, 
O, E, Ne, Na, Mg, and all the heavy neutron-capture elements 
with atomic number in the range [31, 82], including the light-s 
elements (or Is, namely strontium, yttrium and zirconium), the 
heavy-s elements (or hs, namely barium, lanthanum and cerium) 
and lead. The abundance of nitrogen produced in AGB stars is 
uncertain. In AGB stars of mass above about 3Mq, carbon is ef¬ 
ficiently converted to nitrogen at the base of the convective enve¬ 
lope (hot bottom burning, Lattanzio|1991| l. At lower mass, some 
form of deep mixing of the envelope material operates down to 
regions where hydrogen burning occurs, as discussed in Sect.|^ 
but the exact amount depends on the extent of mixing, which is 
very uncertain (e.g.|Hollowell & Iben|1988 [Gallino et'afi 1998[ 


Goriely & Mowlavi|2000[ Stancliffe|2010[ Lugaro et al.|2012 i. 

However, the total amount of Ch-N is conserved, therefore when 
both elements are measured we consider their combined abun¬ 
dances. 

We exclude from Eq. Q the elements from aluminium to 
zinc because they are not involved in AGB nucleosynthesis (alu¬ 
minium is produced only in massive AGB stars), and the dif¬ 
ferences between models and observations reflect a discrepancy 
with our set of initial abundances, as we noticed in the sample 
of C-normal stars. Eor comparison, in some CEMP-s/r stars we 
take into account a smaller set of elements, i.e. only those that are 
mostly produced in AGB nucleosynthesis (C, Mg, Sr, Y, Zr, Ba, 
La, Ce, Pb), as will be explained in Sect. |5.2| Tables [Ai^|A.3| and 
|A.4| in Appendix[A]summarise our results (computed with model 
sets A, B and C, respectively) for the best-fitting models of 43 
CEMP-s stars with number of degrees of freedom v > 2. The 
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number of degrees of freedom is calculated as v = A^obs-4, where 
A^obs is the number of observed elements used in Eq. Q and 4 
is the number of fitted parameters. Our results for the CEMP-s 
stars with v < 1 are shown in Tables |A3|4A.7| 

Some stars fail to match one or more elements within their 
observed uncertainty. To quantify the discrepancy between the 
observations and the model predictions, we calculate the residu¬ 
als, 

7^/ — ^/,obs -'4/,mod ? (2) 




where A,;obs and A,-mod are the observed and modelled number 
abundance of element i of the best fitting model, as in Eq. Q. 

Eor every element i, if our models reproduce correctly the 
nucleosynthesis and mass-transfer process, under the assump¬ 
tion that the measurements are only affected by Gaussian errors, 
the distribution of Ri should resemble a Gaussian function cen¬ 
tred on zero with a standard deviation approximately equal to the 
average observational error. In Eig.[^we plot the distributions of 
Ri computed with the default model set A and binned in intervals 
of width 0.1 for model stars with v > 2. Erom the top left to the 
bottom right we show the distributions of the residuals computed 
for the abundances of carbon (for those stars without a nitrogen 
measurement), Ch-N, sodium, magnesium, strontium, barium, 
europium and lead. The subgroups of CEMP-s/nr, CEMP-s/r 
and CEMP-s/ur stars are indicated with hatched, filled and dot¬ 
ted histograms, respectively. In brackets we indicate the number 
of stars in which the element is measured. 

In Pig. I^we plot, for each star with v > 2, the residuals /?, 
of every observed element. The elements considered in Eq. Q 
are shown as black filled circles, while the other elements are 
shown as grey open circles. The average values of the residuals 
of each element are shown as crosses. Panels (a)-(c) show the 
results obtained with our default model set A for CEMP-s/nr, 
CEMP-s/r and CEMP-s/wr stars, respectively. 

5.1. CEMP-s/nr stars 

Ten of the 67 sampled CEMP-s stars are classified as CEMP- 
s/nr, with [Ba/Ee] > 0.5 and [Eu/Pe] < 1. Eor all the stars 
in this group we find good fits to the abundances. Pour stars 
of this group, CS22942-019, CS22964-161A,B and HD 198269 
belong to binary systems with known orbital periods. In Paper I 
we find the models that best reproduce, at the same time, the 
observed abundances and the orbital periods of these stars. In 
the present work we ignore the constraint on the period in our 
models and we find essentially the same input parameters Mi i, 
M 2 ,i and Mpmz as in the best models of Paper I, whereas there 
are large differences in the initial and final periods. These differ¬ 
ences are due to the amount of mass AM^cc that the secondary 
star has to accrete to reproduce the observed logg and surface 
abundances. Eor example, the model star CS22942-019 accretes 
AMacc ~ 0.3 Mq with all the model sets and consequently the 
of the fit is the same (v^ - 21). On the other hand, the 

initial periods vary by a factor of four (Pi = 6.4 x 10"^ days with 
model sets A and C, Pi = 1.6 x 10“* with model set B), while the 
final periods are Pf = 1.0 x 10^, 3.5 x 10^, Pf = 1.5 x 10“* days 
with model set A, B and C, respectively (observed Porb = 2800 
days). Hence, only model set B reproduces the observed period 
within the uncertainty of our grid of models. 

Eor binary star CS22964-161A,B, which probably formed 
in a triple system (cf. [Thompson et al.|[2008| and Paper I), we 
find the same results as in Paper I without significant differences 
between model sets A, B and C. A binary system of approx¬ 
imately 1.5 Mq accretes a small amount of material (AM^cc ~ 




Fig. 3. Residual distributions in the model-A CEMP-s stars of 
our sample with v > 2. The number of stars in which the element 
is measured is indicated in brackets. The red-hatched, blue-filled 
and grey-dotted histograms indicate the residual distributions for 
CEMP-s/nr, CEMP-i/r and CEMP-i/wr stars, respectively. In 
the top-left panel the residual distribution of carbon is shown for 
those stars without a nitrogen measurement. 


0.06 Mq) from an initially 1.6 Mq primary star in a very wide 
orbit (Pi = 2.7 x 10^, 1.3 x 10^, 2.6 x 10^ days according to 
model sets A, B and C, respectively). Note that the observed 
orbital period in Table A.l does not correspond to the periods 
found in our models, which is the period of the unseen third 
star. The fit to the abundances observed in HD 198269 improves 
compared to the result of Paper I if we ignore the orbital period 
(Porb = 1295 days). With essentially the same assumptions on 
the masses and longer orbital periods (see Tables A.2 A.4 1 our 
model stars accrete AMacc ~ 0.09 Mq and we obtain = 9, 
whereas in our best-fit model in Paper I (with model set B) we 


find AMacc = 0.13 Mq and^^;^ = 18. 

Pig.|^ shows that the abundances of almost all the elements 
involved in AGB nucleosynthesis are well reproduced and typi¬ 
cally the mean of the residuals is close to zero. The abundance 
of sodium is well reproduced in three out of four CEMP-s/nr 
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Fig. 4. Residuals in our 43 modelled CEMP-s stars with a number of degrees of freedom v > 2. Model set A is adopted. From top 
to bottom: CEMP-s/nr, CEMP-s/r and CEMP-s/ur stars. Black filled circles indicate the residuals of the elements included in Eq. 
([^ to calculate grey open circles indicate the other elements, crosses represent the average values of the residuals. 


stars, while it is overestimated by 0.5 dex in CS22964—161B (in 
which cTNa = 0.3 dex). The mean residual of yttrium (atomic 
number 39) is approximately -0.15 dex, lower than for stron¬ 
tium and zirconium (R ^ 0). The difference is smaller than 
the average error in the measurement of yttrium ((Fobs = 0.20 
dex) but indicates that the abundance of yttrium is systemati¬ 
cally overestimated compared to strontium and zirconium. The 
abundances of erbium (atomic number 68) are underestimated 
by 0.6 and 0.2 dex for stars CS22880-074 and HD196944, re¬ 
spectively, while the observational uncertainty is approximately 
0.2 dex in both stars. An offset of about 0.6 dex is also found 
by |Bisterzo et al. (20121 in their analysis of CS22880-074. In 
addition, for this star our best-fitting model and the models of 


[Bisterzo et al.| ( |2012j l well reproduce other elements traditionally 
associated with the r-process, such as europium and dysprosium. 
This evidence possibly suggests that the abundance determina¬ 
tion of erbium in CS22880-074 is affected by observational bias. 

The abundances of lead are underestimated on average by 
0.18 dex, approximately the mean observational uncertainty of 
this element in CEMP-s stars ((Fobs = 0.20 dex). The distribution 
of the residuals in Fig. [^(red-hatched histogram) shows that the 
residuals of lead are always within 2(Fobs- The distributions in 
Fig. show that the abundances of the other elements are typi¬ 
cally reproduced within 2(Fobs- This indicates that on average our 
models of binary evolution and AGB nucleosynthesis reproduce 
reasonably well the abundances of CEMP-s/«r stars, although 
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the results should be interpreted with care because our sample is 
small and we do not constrain the orbital periods of the modelled 
systems. 


5.2. CEMP-s/r stars 


20 stars in our sample are classihed as CEMP-i/r stars. Two 
of these stars, BS16080 - 175 and BS17436 - 058, are not in 
Table |A.21 because only 5 elements are observed and therefore 
V = 1. Four CEMP-s/r stars (CS22948-027, CS29497-030, 
HD224959 and LP625^4) are analysed in Paper I, and we hnd 
that to reproduce the large enhancements observed in most of 
the neutron-capture elements the secondary star has to accrete 
a large amount of mass, AMacc > 0.2 Mq. At the observed or¬ 
bital periods only model set B predicts efficient mass accretion, 
whereas without the constraint on the period also model sets 
A and C predict large mass accretion in much wider orbits, as 
shown in Tables fA^ and A.4 Our best-htting models have rather 
high despite the large accretion efficiency, partly because 
we fail to reproduce the enhanced abundances of the r-elements, 
and this points to a limitation in our nucleosynthesis model, as 
noted also in Paper I. 


For several CEMP-s/r stars we hnd a particularly poor ht, 
with large reduced^ 5, whereas at a visual inspection 
models appear to ht the observed abundances well if Amin/'' - 
To verify if these results are caused only by the discrepancy in 
the r-elements we perform the analysis taking into account 
only nine elements that are produced in large amounts by our 
AGB nucleosynthesis models, and that are frequently detected 
in our observed sample: C, Mg, Sr, Y, Zr, Ba, La, Ce, Pb. This 
choice generally leads to a better ht of the abundances of these 
elements without signihcantly altering the results for the other 
elements. However, even with this choice, in most model stars 
the reduced still exceeds 3, with the exceptions of CS22948- 
027 and CS29497-030. The best-htting models of these two 
stars, calculated taking into account only nine elements, have the 
same primary and PMZ masses as in the models determined in¬ 
cluding all the observed elements (for both stars, Mii = 1.5 Mq 
and MpMz = 2 x 10“^ Mq, with model set A). The secondary 
stars are slightly more massive (by approximately 0.05 Mq for 
both stars with all model sets). In the other CEMP-i /r stars large 
A^;„ are determined because in most of our models there is also 
an issue in reproducing the observed element-to-element ratios. 
In particular, it is difficult to reconcile the large enhancements 
of heavy-s elements and lead with the relatively small abun¬ 
dances of lighter elements (carbon, sodium, magnesium, light-i 
elements). As an example, we compare the modelled and ob¬ 
served abundances in star CS22898-027 in Fig. A model of 
a 0.5 Mq secondary star that accretes AMacc ~ 0.3 Mq from a 
2 Mq primary star with Mpmz = 4 x 10“^ Mq (dotted line) repro¬ 
duces the large enhancement of the heavy-s elements and lead, 
but overestimates carbon by 0.8 dex, sodium and magnesium 
by more than 1 dex and the light-i elements by at least 0.5 dex 
(consequently, A^i„ = 632). The best compromise between the 
light and heavy elements is the model with M] i = 1.5 Mq and 
4 ^pmz = 6.6 X 10 Mq, shown as a solid line in Fig. where 
all elements up to zirconium are overestimated, the neutron- 
rich elements between barium and lead are underestimated and 
A^in = 100 (v = 11). Similar results are found with all model 
sets. 


In Fig. 1^ we plot the hs-to-ls ratio versus the abundance 
of europium observed in our sample CEMP-s/nr and CEMP- 
s/r stars. The abundances of hs and Is are dehned, respectively. 



0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 


atomic number 

Fig. 5. CEMP-s/r star CS22898-027. Points with error bars: ob¬ 
served abundances, in black the elements adopted in Eq. ([ij, in 
grey the other elements. Red solid line: the best-fitting model 
found with model set A, M\ { - 1.5 Mq and Mpmz = 6.66 x 
IO^^'Mq. Blue dashed line: alternative model with Mi,i = 2Mq 
and Mpmz = 4 x 10“^ Mq. Lower panel: the residuals of the two 
models, computed with Eq. Q, are shown as red plus signs with 
error bars and blue triangles, respectively. 


as [hs/Fe] = ([Ba/Fe] -H [La/Fe] -H [Ce/Fe])/3 and [Is/Fe] = 
([Sr/Fe] + [Y/Fe] + [Zr/Fe])/3. If one of these abundances is 
not available we average among the elements present in our 
database. Orange circles and black crosses indicate, respectively, 
the observed stars for which we hnd a ht with small reduced a^ 
(a^ /v < 3) and those for which our best ht has /v > 3. 
The dotted line indicates the maximum value of the hs-to-ls ra¬ 
tio predicted in any of our models, [hs/lslmax = 0.88. We notice 
that; 


[hs/ls] is generally larger for increasingly large [Eu/Fe], be¬ 
cause the abundances of heavy-s elements and europium cor¬ 
relate in metal-poor stars (e.g. Jonsell et al.||2006] l, whereas 
the abundances of light-i elements are essentially indepen¬ 
dent of europium. 

11 stars have observed [hs/ls] < [hs/ls]max. and for all we 
hnd a model with a^ /v < 3. Only 2 out of 11 stars are 


CEMP-s/r stars (CS22881-036 and HD187861), and these 
have relatively low europium enhancements, [Eu/Fe] < 1.4. 

- 17 stars have observed [hs/ls] > [hs/ls]njax; for only two of 
these stars we hnd a model with y^. /v < 3 and one of these 

^min' 

is a CEMP-s/nr star (HD198269). 

- 16 out of 18 CEMP-s/r stars have [hs/ls] > [hs/ls]niax 
(only the CEMP-s/r stars with v > 2 are shown in Fig. 
1^. The value [hs/ls]niax is found in a model AGB star of 
mass Ml i 1.5 Mq: consequently, an initial primary mass 
Mia = 1-5 Mq is selected in almost all our model CEMP-s/r 
stars. However, only three model CEMP-s/r stars reproduce 
the observed abundances well (y^ /v < 3). 


The distributions of the residuals of CEMP-s/r stars in 
Fig. rehect the difficulty in predicting the correct element- 
to-element ratios in most CEMP-s/r stars. The abundances of 
carbon, sodium, strontium, yttrium, and zirconium are on aver¬ 
age overestimated. The mean residual of carbon 'kR - -0.3 dex, 
whereas the mean residual of Ch-N is only R — -0.1 dex. This 
may be interpreted as the effect of extra mixing of protons in 
AGB stars converting carbon to nitrogen. However, the residuals 
of C-hN have a large dispersion, therefore it is not possible to 
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Fig. 6. [hs/ls] vs [Eu/Fe] in our sample CEMP-i/nr (those with 
[Eu/Ee] < 1) and CEMP-i/r stars. Orange hlled circles indi¬ 
cate the stars for which we find a good fit < 3). Black 

crosses indicate the stars poorly reproduced > 3). The 

dotted line indicates the maximum hs-to-ls ratio in our models, 
[hs/ls]tnax = 0 . 88 . 


derive strong conclusions about the efficiency of the mixing pro¬ 
cess. Oxygen is underestimated on average by 0.6 dex, approx¬ 
imately the same as in C-normal metal-poor stars. Because the 
abundance of oxygen is not much affected by AGB nucleosyn¬ 
thesis, this result suggests that the discrepancy could be reduced 
by adopting a larger initial abundance of oxygen in our models. 

The abundance of sodium is always overestimated, on aver¬ 
age by 0.5 dex. In stars HE0338-3945 and HE1305 -h 0007 the 
residuals of sodium are approximately -1 while the observa¬ 
tional uncertainties are 0.1 and 0.2 dex, respectively. The abun¬ 
dance of sodium is increasingly larger the larger the mass of the 
PMZ. In our model CEMP-s/r stars we mostly assume Mpmz ^ 
2 X 10“^ Mq to reproduce the observed abundances of heavy-i 
elements and lead; consequently we overestimate sodium. An 
offset of approximately 0.8 dex in the sodium abundance of 
HE1305-I-0007 is also found by Bisterzo et al. (2012 1 with mod¬ 
els that reproduce the large enhancements of heavy-s elements 
and lead. In contrast, the best ht of |Bisterzo et al.| ( |2011) 1 matches 
the sodium abundance in HE0338-3945 within the observational 
uncertainty. This is probably because the abundances of all r- 
process isotopes are initially pre-enriched, [r/Fe]'"' = 2 dex, 
and no dilution of the accreted material on the secondary star 
is considered. As a consequence, even a low-mass AGB model 
= 1.3 Mq) predicts large abundances of heavy elements, 
while it does not produce much sodium. 

Magnesium is the only light element with mean residual ap¬ 
proximately zero in our models. As shown in Fig. the max¬ 
imum discrepancy between the observed abundances and our 
models is 0.4 dex, within two times the average observational 
uncertainty (cTobs = 0.2 dex). 

The computed residuals of Is, hs and lead conhrm the re¬ 
sults obtained for the example star CS22898-027: the elements 
of the hrst s-peak are overestimated on average by 0.2 - 0.5 dex, 
whereas lead and the elements of the second s-peak are under¬ 


estimated by 0.1 - 0.2 dex, with the exception of cerium and 
lanthanum (R 0). A systematic discrepancy, negative for stron¬ 
tium and positive for barium and lead, is present also in the dis¬ 
tributions of the residuals in Fig. (blue histograms). Because 
none of the model stars reproduces the large observed hs-to-ls 
ratio, to minimise our model predicts larger abundances of 
strontium, and smaller abundances of barium and lead compared 
to the observations. 

The yttrium residuals are typically lower than strontium and 
zirconium, analogously to the results obtained for the sample of 
CEMP-i/nr stars. The abundance of niobium (atomic number 
41) is observed only in star CS29497-030 and is underestimated 
by 0.34 dex (with observational uncertainty 0.2 dex). This is pe¬ 
culiar because according to the models niobium is formed from 
the radioactive decay of ®^Zr, and the abundance of zirconium is 
well reproduced in this stars (R =» 0). This discrepancy is possi¬ 
bly related to the uncertain neutron-capture cross section of ®^Zr, 
as we discuss in Sect. 17. II 

The r-elements are systematically underproduced, perhaps 
not surprisingly because in our nucleosynthesis model the r- 
process is not included. Our models typically produce [Ba/Eu] 
close to 1 , and the maximum europium enhancement is approx¬ 
imately [Eu/Fe] = 1.5, whereas in most CEMP-i/r stars the 
observed [Ba/Eu] is below 0.6 and [Eu/Fe] > 1.5. 


5.3. CEMP-s/ur stars 


37 stars in our sample do not have an observed abundance of 
europium, and therefore are classihed as CEMP-i/wr stars, 15 
of which are listed in Tables IA.211A.4l because at least 6 ele¬ 
ments have been observed. Three of these 15 stars have mea¬ 
sured orbital periods and are discussed in Paper I. One of them, 
BD-i-04°2466, has a very wide orbit (Porb ~ 4600 days), and 
even if we do not consider the period constraint we hnd the 
same input parameters as in Paper I. On the contrary, to re¬ 
produce the period of HD201626 (407 days), our model bi¬ 
nary stars experience a common-envelope phase which shrinks 
the orbit. Without the period constraint the model progenitor 
system of HD201626 has a 1.4 Mq primary star that transfers 
AMacc ~ O.I 8 M 0 to its companion star in a wide orbit (cf. 
Tables A.2 - A.4| l. Consequently, the systems do not experience a 
common-envelope phase and therefore the hnal periods derived 
in the models are 20 to 600 times longer than that observed. 
Without the period constraint we hnd x^ ~ ^ (v = 5), whereas 
in our best-ht model in Paper I we hnd x^ ~ 14 (v = 6 ). In 
Paper I we hnd that the model for star HE0024-2523 needs to 
experience inefficient common-envelope ejection to reproduce 
the large enhancements of the heavy-s elements. If we relax the 
period constraint the system does not enter a common-envelope 
phase and the hnal period of the system is 475 days (with model 
set A), two orders of magnitude longer than the observed period 
of 3.14 days. 

Fig. 1^ shows that carbon, magnesium, zirconium and the 
heavy-s elements are reproduced within 0.1 dex on average. The 
elements with the largest discrepancy with the observations are 
Ch-N (R - -0.45), sodium (R - -0.25), strontium and yt¬ 
trium (R - -0.2) and lead (R - 0.25). The abundance of nitro¬ 
gen is observed only in hve systems (BD-i-04°2466, HD 13826, 
HD187216, HD201626, HD5223) and the average results are 
greatly affected by the poor ht of star HD187216, which is pos¬ 
sibly biased by large observational errors ( [Kipper & Jorgens'm 
19941, and in which our model cannot reconcile the large en¬ 
hancements of the s-elements with the solar ratios observed for 
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Fig. 7. Initial parameters of our 23 modelled stars With < 3 and v > 2. CEMP-s/nr, CEMP-s/r and CEMP-s/ur stars are 

shown respectively with hlled circles, crosses and open diamonds. The bars indicate the conhdence intervals. Prom left to right: 
primary mass, Mi^ mass of the PMZ, Mpmz. secondary mass, M 2 ,i, and orbital period, P^. 


nitrogen, sodium and magnesium. Our best ht reproduces the 
i-elements but overestimates C+N and sodium by 1.1 and 0.7 
dex, respectively. The average residuals of C+N and sodium in 
CEMP-i/i/r stars decrease if we remove HD187216 from the 
sample {R - -0.24 and 0.18, respectively). The distribution of 
strontium is greatly affected by the poor fit of star HE0212- 
0557, where the observed heavy-s elements are enhanced by 2 
dex whereas [Sr/Pe] ^ 0, that is 1.2 dex lower than predicted by 
our best model. Our model also overestimates the abundance of 
yttrium ([Y/Pe] = 0.7) by 0.6 dex, and the abundances of carbon 
and magnesium by 0.5 dex. The models of Bisterzo et al. ( 2012| l 
have similar discrepancies for these elements. If we do not con¬ 
sider this star, the mean residuals of strontium and yttrium are 
-0.1 and -0.15, respectively. Although the differences are small, 
it may indicate the same systematic effect found in CEMP-i/nr 
and CEMP-i/r stars. Pig. shows that the abundances of bar¬ 
ium and lead are reproduced within 0.3 dex, that is 1.5 times the 
average observational uncertainty (cFobs ~ 0.2), with the excep¬ 
tion of the model for star HE0024-2523 that underestimates the 
observed [Pb/Pe] = 3.2 by almost 1 dex. 


6. Confidence limits on the initials parameters 

The results presented in Sect.j^are based on the best-fitting mod¬ 
els found from our ;^f^-minimization procedure. The conhdence 
intervals of the input parameters of our models (i.e. the initial 
primary and secondary masses, the initial period and the PMZ 
mass) are determined with the same procedure as described in 
Paper I. In biref, on our grid of binary models and associated 
we hx one of the input parameters p and for each grid value of p 
we determine the minimum with respect to the other parame¬ 
ters, dehning a function^^(p). A conhdence region is an interval 
of p for which - X^ip) ~ X^m below a certain threshold, 
where is the^^ of the best ht. If the observational errors are 
Gaussian then the thresholds Ax^ = 1,4 and 9 correspond to the 
conhdence intervals with, respectively, 68.3%, 95.4% and 99.7% 
probability that the actual p is in this interval. As we mention in 
Paper I, the observational errors in our sample are not necessarily 
Gaussian, and they may in some cases be affected by systematic 
effects. Consequently, the thresholds Ax^ - 1, 4, 9 should not 
be used to calculate the theoretical Gaussian probabilities, and 


provide only an indication of the conhdence levels. At a visual 
inspection models with x^ below the threshold A)^ - 4 predict 
very similar abundances to the best-htting models. On the other 
hand, models with larger x^ are clearly distinct from the best 
model and the observations are not well reproduced. 

Pig. 1^ shows the input parameters of the 23 model CEMP- 
s stars with /v < 3 and v > 2 as determined with model 
set A of a spherically symmetric wind with WRLOP wind- 
accretion efhciency. CEMP-s/nr, CEMP-s/r and CEMP-s/wr 
stars are indicated as hlled circles, crosses and open diamonds, 
respectively. The horizontal bars represent the conhdence inter¬ 
vals determined with the threshold Ax^ < 4. Model stars with 
X^iJv > 3 are not included because they do not reproduce the 
observed abundances well enough, while stars with v < 2 do not 
provide strong constraints on the models. The conhdence inter¬ 
vals shown in Eig.jTjare summarised in Table A.8 In some stars 
we hnd multiple local minima in the x^ distribution. In these 
stars different combinations of initial parameters result in mod¬ 
els with similar surface abundances and hence almost equal 
As an example, for star CS22880-074 two local minima of x^ 
are found for models with primary mass between 0.9 and 1.1 Mq 
and initial period between 600 and 1.6 x 10"^ days, and another 
around Mi i = 1.5Mo and - [1.6, 3.5] x 10^ days. 

Most of the best-htting models have primary mass below 
2Mq. AGB models with masses in this range produce a signih- 
cant amount of s-process elements regardless of the mass of the 
PMZ because of proton-ingestion events that occur in the hrst 
thermal pulses (as described by Lugaro et al.|2012' which they 
refer to as regime 4 of neutron capture). As a consequence, the 
mass of the PMZ in our model stars is not well constrained, as in¬ 
dicated by the large conhdence intervals in Pig. |7]r. Most model 
stars have Mpmz larger than 10 “^ Mq because in this range low- 
mass AGB models (Mpi < 3 Mq) produce positive [hs/ls] and 
[Pb/hs] as observed in the majority of our sample CEMP-i stars. 

Stars of mass around 1.5 Mq and with Mpmz > 2 x 10~^ Mq 
produce the largest [hs/ls] value and [Pb/hs] > 0. As a con¬ 
sequence, most observed stars with very enhanced heavy-i ele¬ 
ments and lead are best modelled with primary stars of masses 
between Mi; = 1.4 Mq and Mi i = 1.6 Mq. CEMP-s/r stars 
are mostly observed with [hs/ls] close to one or higher, as dis¬ 
cussed in Sect. |5.2[ and therefore our best hts are found with 
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Fig. 8. Initial and final periods, Pi and Pf (respectively, top and bottom panels), in our models with;^;'^ /v < 3 and v > 2. Symbols 
are the same as Fig. [7] From left to right model sets A, B and C are shown (panel a is the same as Fig?j^). 


M\;i - 1.5 Mq. Two CEMP-i/nr stars are modelled with 
in the mass range [1.75, 3] Mq, in which the [hs/ls] ratio varies 
approximately between -0.5 and 0.5 for MpMz betwee n 0 and 
2x 10“^ Mq while [Pb/hs] is always positive (regime 2 of|Lugaro 


et al. 2012|l. The secondary mass (Fig. [7};) in the models de¬ 


pends on the observed enhancements. The default assumption in 
our models, as in Paper I, is that accreted material is efficiently 
mixed by thermohaline mixing. As a consequence, a large de¬ 
gree of accretion is required when the observed abundances are 
greatly enhanced. Hence, the secondary star has to be initially 
less massive to reach a mass of approximately 0.8 - 0.9 Mq after 
accretion, and still be visible at f > 10 Gyr. 


The input parameters Mi i, MpMz M 2 ,i are similar in all 
model sets. On the contrary, there are several differences be¬ 
tween the initial and final periods, as shown in the upper and 
lower panels of Fig. In model set C the binary systems typi¬ 
cally start in wider orbits compared to the other model sets, be¬ 
cause the dependence of the WRLOF process on the separation 
favours wider orbits compared to model set B and because of the 
more efficient angular momentum loss compared to model set A. 
The hnal periods computed with model set A and C are similar 


in many stars because the mass-transfer algorithm is the same 
in the two models, and therefore the binary stars need to have 
approximately the same separation during accretion to transfer 
the same amount of material. The majority of stars in model 
sets A and C have final periods longer than 4,600 days, which 
is approximately the longest period measured in our sample. On 
the contrary, in model set B mass transfer is more efficient in 
close orbits, and consequently the majority of the systems are 
predicted at periods shorter than 10,000 days. This is the only 
model set that predicts signihcant mass accretion in close or¬ 
bits, as shown also in Paper I, and ten modelled stars have fi¬ 
nal periods less than 4,600 days. The stars that need to accrete 
a small amount of mass to reproduce the observed abundances 
have longer periods because of the enhanced wind-accretion rate 
the and efficient angular momentum loss. 
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7. Discussion 

7.1. Comparison between observed and modelled 
abundances 

Our best-fitting models are found with initial primary mass in 
the range between Q.9Mq and 3Mq. In this mass range the neu¬ 
tron source is the *^C(a, n)'®0 reaction that operates on pro¬ 
duced by the inclusion of a PMZ or by ingestion of protons in 
the stellar interior during the thermal pulses (regimes 2 - 4 in 
Lugaro et al.|2012]l. The abundances of ^-elements predicted in 


the models depend on the mass of the AGB star and of its PMZ. 
We note that: (/) [Ba/Eu] is always close to 1 dex, (ii) the value 
of [hs/ls] varies from [hs/ls] « -0.5 for Mp^z - OMq up to 
[hs/ls] ai 0.9 for MpMz > 2 x IO^^Mq and Mi i - I.SMq, and 
{Hi) [Pb/hs] is always positive, and in some cases greater than 
one, for Mi j above approximately 2Mq. 

The success of our models in reproducing the observed abun¬ 
dances varies significantly for different classes of stars. CEMP- 
sjnr stars typically exhibit [hs/ls] 0.5 dex and [Ba/Eu] be¬ 
tween 0.5 and 1 dex, therefore our models generally well repro¬ 
duce the abundances of the s-elements and also the r-elements. 
On the other hand, most CEMP-s/r stars have [hs/ls] > 1, 
and therefore our models are not able to reproduce the light-i 
peak and the heavy-i peak simultaneously. In most cases the 
model with the minimum systematically overestimates the 
light-s elements, on average by 0.3 dex, and underestimates the 
heavy-s elements, on average by 0.1 dex. Eurthermore, generally 
[Ba/Eu] < 1 and the abundances of all r-elements are underesti¬ 
mated by a factor of two to up to one hundred. In our CEMP-i/wr 
stars the abundances of the heavy-i elements are typically re¬ 
produced within the observational uncertainty, as in CEMP-s/nr 
stars, whereas the light-s elements are systematically overesti¬ 
mated on average by 0.15 dex, similar to CEMP-s/r stars. These 
results are probably related to the reason why the abundance of 
europium is indeterminate. This group is likely to contain a mix¬ 
ture of CEMP-i/nr stars and CEMP-i/r stars. Hence, in some 
stars the abundance of europium is actually low, whereas in other 
stars europium is enhanced but it is not detected, for example be¬ 
cause the europium lines are blended. 

The abundance ratio between carbon and heavy-s elements is 
generally positive in our AGB-nucleosynthesis model, [C/hs] > 
0. Consequently, in our synthetic stars high abundances of 
heavy-i elements are associated with large enhancements of car¬ 
bon. A similar correlation is also found in our observed sam¬ 
ple: compared to CEMP-s/nr stars, CEMP-s/r stars typically 
exhibit higher abundances of carbon (Pig. [T]) and heavy-s ele¬ 
ments (as indicated by the increase of [hs/ls] with [Eu/Pe] in 
Pig. |6]l. However, our models overestimate the carbon-to-hs ra¬ 
tio in CEMP-i/r stars. As a consequence, in the majority of 
these stars the carbon abundance is overestimated by our mod- 


els with high [1 
by Izzard et al. 

is/Pe]. 
([2009 

A similar problem is also encountered 
e.g. Pig. 13b) and in some models of 

Bisterzo et al. ( 

2011[|2012 

1 . This probably indicates that the 
s should be able to produce higher 

AGB-nucleosynthesis mode 


hs abundances for a given amount of carbon, as we also discuss 
in a forthcoming paper (Abate et al.|2015b in press). 


is very uncertain ( |Milam et al.||2009l |Stanchff^|2010[ [Karakas] 
|et al.|20T0| . If we take into account the combined abundance 
of carbon and nitrogen our results improve signihcantly, and the 
observations are always reproduced within at most twice the ob¬ 
servational uncertainty. 

Oxygen is underestimated in almost all stars of our sample, 
on average by 0.4 dex, approximately the same amount as in C- 
normal stars (cf. Sect. |^. The abundance of oxygen does not 
change much during AGB nucleosynthesis, and most likely the 
discrepancy could be solved by adopting a larger initial abun¬ 
dance without signihcantly affecting the abundances of the other 
elements. 

The abundance of sodium is overpredicted in all our model 
stars, and this points to a general issue in our nucleosynthe¬ 
sis model. Sodium is produced in the intershell region of the 
AGB star due to proton and neutron capture on the abundant 
^^Ne. The abundance of sodium increases rapidly with increas¬ 
ing PMZ mass and therefore models with a smaller PMZ pre¬ 
dict lower abundances of sodium at each stellar mass. However, 
also the abundance of magnesium is sensitive to the mass of the 
PMZ and a large PMZ is required in most stars to reproduce 
its abundance, as well as the abundance of heavy-i elements. 
Lugaro et al.|(201^ suggest three effects that may help lower 


the predicted sodium abundance: (i) in the detailed models of 
AGB nucleosynthesis the density prohle of protons introduced 
to m ake the ^^C-pocket could de crease more rapidly with depth 
(e.g., Goriely & Mowla^|2000| , (ii) a smaller ^^Ne(p, yj^^Na 
reaction rate or (Hi) a larger ^^Ne(/7, a)^°Ne reaction rate, both 
of which are very uncertain ( [Iliadis et al.|2010| l. 

In many CEMP-i stars the abundance of yttrium is system¬ 
atically lower than strontium and zirconium by 0.1 - 0.2 dex, 
whereas in our models these elements are generally produced 
in approximately the same amount for every stellar mass and 
PMZ. This suggest that less yttrium should be produced in the 
models. Yttrium is produced by the reaction **Sr(n, yj^^Sr and 
the subsequent decay of ®®Sr to with a half life of 51 days. 
The neutron-capture cross section of *^Sr is uncertain by a fac¬ 
tor of two ( |Bao et al.|2000 1 . Qualitatively, a larger cross section 
implies that before the decay of ^®Sr the reaction ^®Sr(n, y)®°Sr 
may occur and after two consecutive /1-decays ®°Zr is produced 
while *®Y is skipped and this effect may help to reduce the dis¬ 
crepancy. 

Niobium is observed in the stars CS29497-030 and 
HD187216 and it is underestimated by 0.35 dex and 0.69 dex, 
respective ly (although the abundanc es of HD187216 are very 


uncertain. Kipper & Jorgensen] 1994|l. However, ®^Nb is formed 


Some discrepancies between models and observations are 
found in all three classes of CEMP-s stars. The abundance of 
nitrogen is underestimated in 15 out of 16 stars by 0.1 to 1.8 
dex. The difficulty in reproducing the observed abundance of ni¬ 
trogen is a well known issue of AGB nucleosynthesis models. 
In AGB stars of mass below approximately 3Mo some nitrogen 
is produced from carbon by a deep mixing process that operates 
at the bottom of the convective envelope, but the exact amount 


by jS-decay of ^■’Zr and in the two stars zirconium is well repro¬ 
duced. A similar problem is found by |Kashiv et al. (2010 1 in the 
context of presolar SiC grains, in which the Nb/Zr ratio is sys¬ 
tematically higher than in Cl chondrites. These authors suggest 
that a cross section of ®^Zr reduced by a factor of two removes 
the discrepancy and we note that this solution would help to re¬ 
duce the difference observed in the two CEMP-s stars. New mea¬ 
surements of the ®^Zr cross section are currently being analysed 
and will provide a more accurate determination of the reaction 
rate ( Lugaro et al.|2014| . 

The abundances of barium are underestimated by our models 
in all but two CEMP-s/r stars and in the majority of the CEMP- 
s/wr stars. However, we note that the observed barium abun¬ 
dances might be affected by large spectroscopic uncertainties 
because most of the observed barium lines are strong resonance 
lines, often saturated and sensitive to non-LTE effects ( |Busso| 
|et al.| 1995) [Andrievsky et al.|2009[|Masseron et al.|2010|l. 
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The observed abundance of lead in most CEMP-i stars is 
larger than predicted by the models. The lead abundance in¬ 
creases with increasing stellar mass ([Pb/Fe] is maximum for 
M\ i between 2Mq and 3Mo) and PMZ mass. In most of the stars 
in our sample the choice of a more massive primary star im¬ 
plies that the abundances of sodium, magnesium, and the light-i 
elements are overestimated. This suggest that nucleosynthesis 
models of low-mass stars (Mi,i < 2Mq) should produce higher 
[Pb/Na], [Pb/Mg] and [Pb/ls]. On the other hand, the determina¬ 
tion of lead abundances generally rely on one single line blended 
by CH-lines, therefore the observed lead abundances may be af¬ 
fected by systematic errors ( Masseron et al.|2010| l. 

Lead is one of the most neutron-rich elements produced by 
AGB stars, and therefore its abundance is strongly metallicity- 
dependent ( [Travaglio et al.|20OT I. To assess the effect of metal- 
licity on our results, we use the theoretical models presented in 
Fig. 1 of |Travaglio et al.| ( |2001| l to construct a table of correction 
factors which account for the differences in lead abundance pre¬ 
dicted by their models in the range —2.8 < [Fe/H] < —1.8 and 
our model at [Fe/H] = -2.3. We correct the lead abundances 
predicted in our models by interpolating in this table the iron 
abundances observed in CEMP-s stars. However, the distribu¬ 
tion of the residuals of [Pb/Fe] does not improve with this anal¬ 
ysis, even if we vary the correction factors based on different 
models of Travaglio et al. ( 2001| l. The residuals of [Pb/Fe] com¬ 
puted with our models do not exhibit a dependence on [Fe/H]. 
Consequently, even correcting for different metallicities, the dis¬ 
crepancies between our model predictions and the observations 
are not reduced. Hence, as yet it is not possible to quantify the 
errors introduced in our analysis by the assumption of models 
with fixed metallicity. A larger sample of reliable abundances of 
lead (and other i-elements) in CEMP-i stars will help to address 
this issue. 


7.2. Comparison with previous resuits 


We compare our results with the analysis performed by Bisterzo 
et al. (20121 on 94 CEMP-s stars in the range of metallicity 


-3.5 < [Fe/H] < -1.0, 60 of which have -2.8 < [Fe/H] < 
-1.8 and are in common with our sample. One CEMP-i/r star 
(HE2258-6358, analysed by |Placco etal.|2013| ) and six CEMP- 
sjur stars are only included in our sample. In five of these stars 
the abundances of only a few elements have been determined 
(hence v < 1), whereas HD 187216 was not included in the study 
of Bisterzo et al.| ( 2012| l because of the possibly large observa¬ 
tional biases. 

Bisterzo et al. p012|l compare the observations with the 


results of AGB nucleosynthesis models of single stars with 
masses between 1.2Mq and 2Mq in a range of metallicities 
and '^C-pockets. A quantitative comparison with their results 
is not straighforward because of the intrinsic differences in the 
methods used in the two studies. For example, the treatment 
of the '^C-pocket/PMZ is not the same in the AGB models 
by Bisterzo et al. ( 2012| l and by Karakas (2010 on which our 
model is based). In the former, the parameter that determines 
the efficiency of the s-process is the abundance of in the 
pocket, which is added at the top of the intershell region at 
each third dredge-up. In the latter, the size of the PMZ deter¬ 
mines the mass of a layer of free protons that are mixed with 
the intershell. The protons react with the elements in the in¬ 
tershell and the mass fractions of '^C, and other isotopes 
are calculated consistently. The standard mass of the '^C-pocket 
in the models by Bisterzo et al. \2Q\2\ corresponds roughly to 
MpMz - 2x IO^^Mq in our models. Another important differ¬ 


ence of our study is that we explicitly compute the binary evo¬ 
lution, wind mass-transfer process and the mixing effects while 
in the paper by |Bisterzo et al. ( |2012| l these aspects are mimicked 
by including a dilution factor, i.e. the ratio between the amount 
of accreted material and the envelope mass of the accreting star. 

Despite these differences, in 46 of the 60 CEMP-s stars com¬ 
mon to both studies we find primary masses that are consistent 
with the results of [Bisterzo et al.|P0121 l within the confidence 
intervals of the models. Among these stars are 21 of the 23 stars 
for which we find a good fit (with ■ Iv < 3 and v > 2). The 
two exceptions are HD5223 and HE2232-0603. Compared to 
the results of [Bisterzo et al.j ( |2012| ), in these CEMP-s/ur stars 
our models with low primary mass (0.9 - 1.1 Mq) better match 
the observed carbon and lead (in both stars), fluorine, sodium 
and magnesium (in HD5223), whereas larger ^ overestimate 
the observed abundances. 

Of the stars for which there is no agreement with the masses 
predicted by Bisterzo et al. ( |2012 1 , ten are CEMP-s/r stars, 
while the other four are CYMP-sjur stars. One possible expla¬ 
nation of these discrepancies is in the way [Bisterzo et al.] ( [2012[ l 
account for the contribution of the r-process in these stars. The 
authors assume that CEMP-i/r stars were formed from gas pol¬ 
luted by the explosion of a Type II supernova, and therefore 
adopt in their models an initial enrichment of the r-process com¬ 
ponent of the elements heavier than iron, between [r/Fe]“‘ = 0.5 
and [r/Fe]™ = 2. In CEMP-s/r stars this initial enhancement is 
chosen to reproduce the observed [Eu/Fe]. The enhanced abun¬ 
dance of the r-process component does not alter significantly the 
abundance of the i-process elements, except in low-mass stars 
(Ml i < 1.7Mo) with [r/Fe]‘“ = 2, in which case the heavy-s 
elements are enriched by an extra factor of a few tenths of a dex 
(e.g.. Fig. 6 of [Lugaro et al.[[2012|l. As a consequence, the r- 
enriched models of [Bisterzo et al.[ ( [2012[) with low initial masses 
(M = 1.2-1.4Mo) and relatively smalP^C-pockets, correspond¬ 
ing roughly to Mpmz ^ IO^^Mq, reproduce the abundance of the 
heavy-s elements, r-elements and lead without greatly overesti¬ 
mating the abundances of the light-s elements. In contrast, in our 
models neutron-capture elements are entirely produced by the 
i-process during the AGB phase. Consequently, our models of 
CEMP-i/r stars are constrained towards relatively massive pri¬ 
mary stars, 1.5Mo < { < 2.9 Mq, with Mpmz > 2 x 10“^ Mq, 

which predict large abundances of elements heavier than barium 
but typically overestimate the observed abundances of lighter el¬ 
ements. 

The origin of the r-elements in CEMP-s/r stars is unclear. 
Many hypotheses have been suggested that can be classified in 
three main classes. In the first class of scenarios the binary sys¬ 
tems were formed in a molecular cloud already enriched in r- 
elements by the nearby explosion of one or more Type II super¬ 


novae ( Cohen et al.|2003[ Jonsell et al.|2006[[Sneden et al.|2008] 


[Bisterzo et al.|20lT 2012[ ). In the second class of scenarios, the 

primary star of the binary system is relatively massive and un¬ 
dergoes the AGB phase, producing the i-elements, and then ex¬ 
plodes as an electron-capture or a Type 1.5 supernova, providing 
the r-elements ( |Zijlstra|2004[|Wanajo et al.|2006[ [Jonsell et al.[ 


20061. If the s-process and the r-process enrichments are inde¬ 


pendent, as suggested in the first class of scenarios, it is difficult 
to explain the correlation between the abundances of barium and 
europium that is observed in CEMP-s/nr stars and in CEMP- 
s/r stars. Moreover, the massive AGB stars invoked in the sec¬ 
ond class of scenarios are rare in a solar-neighbourhood initial 
mass function ( Pols et al.|2012 1 . Also, massive AGB stars pro¬ 
duce more nitrogen than carbon (e.g. Karakas & Lattanzio|2007 1 
and would therefore produce nitrogen-enhanced metal-poor stars 
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rather than CEMP stars. For a thorough discussion of these sce¬ 
narios we refer to|Jonsell et ah (2006 |l and|Lugaro et al. (|2009 1 . 


Lugaro et al. (20091 suggest a third speculative scenario in which 


r- and s-elements are both formed in low-mass extremely metal- 
poor AGB stars ([Fe/H] < -3), that may be able to produce very 
large densities of free neutrons when protons are ingested in the 


region of the He-flash (|Campbell||2007 

Campbell & Lattanzio| 

2008 1 Herwig et al.|201 1| Reifarth et al. 

2014jl. 


7.3. Constraints on binary mass transfer 


We considered several model sets with different combinations 
of assumptions for the models of wind mass-accretion rate and 
angular momentum loss. Because in this study we do not con¬ 
sider constraints on the orbital period there are no significant 
differences between the results obtained with the three model 
sets, with the exception of the distributions of initial and hnal 
orbital periods of the binary systems. Without the constraint on 
the orbital period our fitting procedure mostly selects model sys¬ 
tems with periods longer than 5,000 days, which is more than 
the longest period currently observed. In our default model set 
A relatively long periods are necessary to accrete large amounts 
of mass. Therefore, when large mass accretion is required to 
reproduce the observed abundances we find modelled CEMP- 
s stars with periods approximately 20 - 500 times longer than 
predicted in Paper I, in which the observed orbital periods are 
used to constrain our models (e.g. the stars CS22942 - 019 and 
CS29497 - 030). Shorter periods are found in a few cases if 
the observed abundances are reproduced by a model with a low- 
mass primary star (Mi_i < I.IMq) and small accreted mass. For 
these stars we compute orbital periods of a few thousand days, 
consistent with the results of Paper I (e.g. the stars BD-i-04°2466 
and HD198269). 

With our model set A we hnd that our best-htting model sys¬ 
tems have orbital periods between 600 days and 400,000 days 
with an average period of 90,000 days. With model sets B and 
C we obtain on average Pf - 17,000 and Pf — 60,000 days, 
respectively. On the other hand, the average orbital period of 
the systems analysed in Paper I is approximately 1,500 days. 
Furthermore, a recent statistical study of currently available ra¬ 
dial velocity variations in CEMP-i stars indicates that these stars 
are expected to have a maximum period of about 10,000 days 


and an average period of 400 days (Starkenburg et al. 20141. This 
suggests that a large proportion of the CEMP-s stars in our sam¬ 
ple have periods of less than a few thousand days, i.e. in closer 
orbits than we hnd in our best hts. If this hypothesis is correct, 
then AGB stars in binary systems have to transfer wind mate¬ 
rial with great efficiency at relatively short separations and this 
process is currently possible in our models only with some ad 
hoc assumptions, namely highly efficient mass-transfer and ef¬ 
ficient angular momentum loss. This hypothesis is supported by 
the period distribution derived in barium stars, the equivalent of 
CEMP-i stars at solar metallicity (e.g. Jorissen et al.|1998| l, and 
by observations of ongoing mass transfer in post-AGB binaries 
(e.g. [Gorlova et al.|2012] l. 


8. Summary and conclusions 

In this work we have used our model of binary evolution and nu¬ 
cleosynthesis to calculate the best-fitting models to the surface 
abundances observed in a sample of 67 CEMP-i stars. We dis¬ 
tinguish three classes of CEMP-s stars based on the europium 
abundance. In CEMP-s/nr stars the abundance of europium rel¬ 
ative to iron is up to ten times as high as in the Sun. In CEMP- 


s/r stars the europium-to-iron ratio is higher than ten times the 
solar ratio. In CEMP-s/ur stars the abundance of europium is 
indeterminate. From the comparison between the observed and 
modelled abundances of these stars we draw the following con¬ 
clusions. 

- The model stars that provide the best fit to the abundances 
observed in CEMP-i stars have low initial mass (up to 
2.5Mo). In this range, the main neutron source for the s- 
process is the *^C(a, n)'®0 reaction. 

- The chemical properties observed in CEMP-s/nr and 
CEMP-s/r stars are fundamentally different. The results of 
our models are consistent with the abundances observed 
in CEMP-i/nr stars. On the contrary, most of the model 
CEMP-i/r stars overestimate the observed abundances of 
light-i elements and underestimate the abundances of heavy- 
s elements, r-elements and lead. In CEMP-s/r stars the ra¬ 
tio of heavy-s elements to light-s elements is too high to be 
reproduced in our models. This result points to a different 
nucleosynthesis process at the origin of the abundances in 
CEMP-i/r stars, that enhances the heaviest elements with¬ 
out affecting the abundances of elements such as sodium, 
magnesium and light-s elements. 

- The abundance of sodium is always overpredicted in our 
models. This discrepancy likely points to a problem in the 
numerical treatment of the partial mixing zone. A proton pro¬ 
file in the partial mixing zone weighted towards a low proton 
abundance may contribute to reducing this discrepancy. 

- The orbital periods predicted in this work are significantly 
longer than the results of Paper I, in which our models are 
constrained to reproduce the observed orbital periods of the 
systems. This indicates that either the currently known pe¬ 
riods of CEMP-i stars are biased towards close orbits and 
are not representative of the whole population, or our model 
of wind accretion in binary stars needs to efficiently transfer 
mass and lose angular momentum in close orbits. 
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Appendix A: Tables 

In this section we summarise the observed and modelled proper¬ 
ties of the 67 CEMP-i stars in our sample. Table |A.l1 shows the 
observed temperatures, surface gravities and abundances of se¬ 
lected elements. The parameters of the best-htting models to the 
CEMP-s stars with at least of two degrees of freedom (v > 2), 
as computed with model sets A, B and C, are shown in Tables 
|A.2HA.4| as follows. 


Columns 2-5: the htted parameters Mij, Mpmz, A/ 2 ,i and P{. 
Columns 6-8: the mass accreted by the secondary star, 
AMacc, the orbital period of the binary when the secondary 
star best reproduces the observed logg and surface abun¬ 
dances, Pf, and for stars with observed orbital periods, the 
hnal period determined in Paper I. 

Columns 9-11: number of degrees of freedom of 

the ht, V, and the reducedof the best ht, i.e.;^f^ /v = 


In Tables A.5 A.7| the physical parameters of the best-htting 
models of stars with v < 1 are shown.The conhdence intervals 
of the input parameters of our models with ■ / v < 3 and v > 2 


computed with model set A are shown in Table A.l 
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Table A.l. Surface gravities, temperatures and chemical properties observed in the 67 CEMP-s stars of our sample. 


ID 

logiote/cms 

Teff/K 

Porb/days 

number of 
observed elements 

[Fe/H] 

[C/Fe] 

[VFe] 

[Eu/Fe] 

CEMP-i/nr stars 

CS22880-074 

3.9 ±0.13 

5917 


21 

-1.9 

1.5 

1.3 

0.6 

CS22942-019 

2.2 ± 0.4 

4967 

2800 

18 

-2.7 

2.2 

1.8 

0.8 

CS22964-161A 

3.7 ± 0.2 

6050 

252.481 

21 

-2.4 

1.6 

1.4 

0.7 

CS22964-161B 

4.1 ±0.4 

5850 

252.481 

22 

-2.4 

1.4 

1.3 

0.7 

CS30301-015 

0.8 ±0.1 

4750 


17 

-2.7 

1.7 

1.5 

0.3 

HD196944 

1.8 ±0.1 

5234 


27 

-2.4 

1.4 

1.2 

0.3 

HD198269 

1.3 ±0.25 

4800 

1295 

11 

-2.2 

1.7 

1.7 

0.9 

HE0202-2204 

1.6 ±0.1 

5280 


21 

-2.0 

1.2 

1.4 

0.5 

HEl 135+0139 

1.8 ±0.1 

5487 


20 

-2.4 

1.1 

1.1 

0.4 

HE2158-0348 

2.5 ±0.1 

5215 


18 

-2.8 

2.1 

1.6 

0.9 


CEMP-i/r stars 


BS16080-175 

3.7 ± 0.2 

6240 


6 

-1.9 

1.8 

1.6 

1.1 

BS17436-058 

2.7 ± 0.2 

5690 


7 

-1.8 

1.6 

1.7 

1.2 

CS22881-036 

4.0 ±0.1 

6200 


14 

-2.1 

2.1 

1.9 

1.0 

CS22898-027 

3.7 ± 0.28 

6110 


22 

-2.3 

2.0 

2.3 

2.0 

CS22948-027 

1.8 ±0.4 

4800 

426.5 

21 

-2.5 

2.4 

2.4 

1.9 

CS29497-030 

4.0 ± 0.5 

6966 

344 

33 

-2.5 

2.4 

2.3 

1.7 

CS29526-110 

3.2 ±0.1 

6500 


18 

-2.4 

2.3 

2.1 

1.8 

CS31062-012 

4.2 ± 0.38 

6099 


24 

-2.8 

2.3 

2.1 

1.6 

CS31062-050 

2.9 ± 0.24 

5489 


37 

-2.5 

1.9 

2.4 

2.0 

HD187861 

2.0 ± 0.35 

4960 


14 

-2.4 

2.0 

1.9 

1.3 

HD224959 

1.9 ±0.25 

5050 

1273 

14 

-2.1 

1.8 

2.2 

1.7 

HE0131-3953 

3.8 ±0.1 

5928 


16 

-2.7 

2.5 

2.2 

1.7 

HE0143-0441 

4.0 ± 0.35 

6305 


22 

-2.4 

2.0 

2.4 

1.7 

HE0338-3945 

4.1 ±0.33 

6161 


32 

-2.5 

2.1 

2.4 

2.0 

HEl 105+0027 

3.5 ±0.1 

6132 


16 

-2.5 

2.0 

2.4 

1.9 

HE1305+0007 

1.5 ±0.5 

4655 


21 

-2.2 

2.1 

2.6 

2.2 

HE2148-1247 

3.9 ±0.1 

6380 


25 

-2.4 

2.0 

2.4 

2.0 

HE2258-6358 

1.6 ±0.1 

4900 


31 

-2.7 

2.4 

2.3 

1.7 

LP625-44 

2.6 ±0.3 

5500 

> 4383 

31 

-2.8 

2.3 

2.8 

1.9 

SDSSJ0912+0216 

4.5 ±0.1 

6500 


28 

-2.6 

2.3 

1.6 

1.3 


CEMP-i/wr stars 


BD+04°2466 

1.8 ±0.2 

5032 

4592.7 

20 

-2.1 

1.3 

1.6 

CS22945-017 

3.8 ±0.2 

6400 


6 

-2.5 

2.3 

0.6 

CS22956-028 

3.9 ±0.1 

6900 

1290 

11 

-2.1 

1.9 

0.4 

CS29509-027 

4.2 ±0.1 

7050 

196 

5 

-2.1 

1.5 

1.3 

CS30338-089 

2.1 ±0.1 

5000 


10 

-2.5 

2.1 

2.3 

HD13826 

0.1 ±0.35 

3540 


9 

-2.5 

1.6 

1.4 

HD187216 

0.4 ± 0.4 

3500 


14 

-2.5 

1.4 

2.3 

HD201626 

2.2 ± 0.25 

5190 

407 

11 

-2.1 

2.1 

2.0 

HD5223 

1.0 ±0.25 

4500 


18 

-2.1 

1.6 

1.9 

HE0012-1441 

3.5 ±0.1 

5730 


7 

-2.7 

1.9 

1.3 

HE0024-2523 

4.3 ±0.1 

6625 

3.14 

17 

-2.7 

2.1 

1.6 

HE0206-1916 

2.7 ±0.1 

5200 


8 

-2.2 

2.1 

2.0 

HE0212-0557 

2.1 ±0.4 

5075 


15 

-2.3 

1.9 

2.2 

HE0231-4016 

3.6 ±0.1 

5972 


17 

-2.1 

1.3 

1.5 

HE0322-1504 

0.8 ±0.1 

4460 


5 

-2.0 

2.4 

2.8 

HE0430-4404 

4.3 ±0.1 

6214 


15 

-2.1 

1.3 

1.6 

HE0441-0652 

1.4 ±0.1 

4900 


11 

-2.6 

1.4 

1.2 

HE0507-1430 

0.8 ±0.1 

4560 

446 

5 

-2.4 

2.7 

1.3 

HE0534-4548 

1.5 ±0.3 

4250 


4 

-1.8 

1.5 

0.6 

HE1045-1434 

1.8 ±0.1 

4950 


5 

-2.5 

3.3 

3.0 

HEl 157-0518 

2.0 ±0.1 

4900 


8 

-2.4 

2.2 

2.2 

HE1429-0551 

1.5 ±0.1 

4700 


7 

-2.6 

2.3 

1.7 

HE1430-1123 

3.8 ±0.1 

5915 


15 

-2.7 

1.8 

1.8 

HE1434-1442 

3.1 ±0.4 

5420 


15 

-2.5 

2.1 

1.3 

HE1443+0113 

1.9 ±0.1 

4945 


5 

-2.2 

2.0 

1.4 

HE1447+0102 

1.7 ±0.1 

5100 


7 

-2.5 

2.5 

2.8 

HE1523-1155 

1.6 ±0.1 

4800 


6 

-2.2 

1.9 

1.8 

HE1528-0409 

1.8 ±0.1 

5000 


6 

-2.7 

2.4 

2.4 

HE2150-0825 

3.7 ±0.1 

5960 


18 

-2.0 

1.4 

1.7 

HE2221-0453 

0.4 ±0.1 

4400 


7 

-2.3 

1.8 

1.8 

HE2227-4044 

3.9 ±0.1 

5811 


14 

-2.4 

1.7 

1.4 

HE2228-0706 

2.6 ±0.1 

5100 


8 

-2.5 

2.3 

2.6 

HE2232-0603 

3.5 ±0.1 

5750 


17 

-2.1 

1.6 

1.6 

HE2240-0412 

4.3 ±0.1 

5852 


13 

-2.2 

1.4 

1.4 

SDSS0924+40 

4.0 ± 0.3 

6200 


12 

-2.6 

2.7 

1.9 

SDSS1707+58 

4.2 ± 0.3 

6700 


8 

-2.6 

2.1 

3.5 

SDSS2047+00 

4.5 ± 0.3 

6600 


12 

-2.1 

2.0 

1.6 


Notes. Barium abundance is used as indicator of s-elements. In HD198269, HD13826 and HD201626 lanthanum is used because barium is not 
observed. 
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Table A.2. Physical parameters of the model CEMP-i stars with number of degrees of freedom v > 2 computed with model set A 
with WRLOF wind-accretion efficiency and spherically symmetric wind. 


ID 


Mfuz 

M2.i 


AMacc 

Pf 

Pf-Paper I 

■ 

xv mm 

V 

X'^ ■ h 

XI min' 

CEMP-j/«r stars 

CS22880-074 

1.1 

1.00 X 10-" 

0.74 

4.54 X 103 

0.09 

6.03 X 103 


16.8 

10 

1.7 

CS22942-019 

1.8 

2.00 X 10-* 

0.54 

6.38 X 10* 

0.30 

1.03 X 103 

2.87 X 103 

21.1 

9 

2.3 

CS22964-161A 

1.6 

2.00 X 10-3 

0.79 

2.74 X 103 

0.04 

5.60 X 103 

5.60 X 103 

19.5 

8 

2.5 

CS22964-161B 

1.6 

2.00 X 10-3 

0.71 

2.74 X 103 

0.02 

5.60 X 103 

5.60 X 103 

11.8 

9 

1.3 

CS30301-015 

2.5 

4.00 X 10-3 

0.76 

4.81 X 103 

0.06 

1.68 X 10* 


11.9 

8 

1.5 

HD 196944 

1.4 

1.00 X 10-3 

0.81 

5.84 X 103 

0.06 

1.14 X 10* 


13.8 

12 

1.2 

HD 198269 

1.0 

4.00 X 10-3 

0.79 

1.63 X 103 

0.09 

1.97 X 103 

1.43 X 103 

9.2 

5 

1.8 

HE0202-2204 

0.9 

3.00 X 10-3 

0.79 

1.19 X 103 

0.06 

1.43 X 103 


5.2 

7 

0.7 

HEl 1354-0139 

0.9 

1.00 X 10-3 

0.79 

3.35 X 103 

0.08 

4.01 X 103 


2.8 

6 

0.5 

HE2158-0348 

1.5 

3.00 X 10-3 

0.69 

1.32 X 103 

0.17 

2.52 X 103 


16.4 

7 

2.3 


CEMP-j/r stars 


CS22881-036 

1.5 

6.66 X 10-* 

0.66 

1.32 X 10= 

0.15 

2.57 X 103 


10.3 

6 

1.7 

CS22898-027 

1.5 

6.66 X 10-* 

0.51 

2.44 X 10* 

0.31 

2.89 X 10* 


99.8 

11 

9.1 

CS22948-027 

1.5 

2.00 X 10-3 

0.64 

9.42 X 10* 

0.24 

1.58 X 103 

4.76 X 103 

29.5 

9 

3.3 

CS29497-030 

1.5 

2.00 X 10-3 

0.46 

2.47 X 10* 

0.28 

2.87 X 10* 

4.04 X 103 

126.0 

16 

7.9 

CS29526-110 

1.5 

2.00 X 10-3 

0.54 

2.42 X 10* 

0.32 

2.92 X 10* 


22.4 

6 

3.7 

CS31062-012 

1.5 

2.00 X 10-3 

0.56 

6.79 X 10* 

0.25 

1.00 X 103 


17.5 

8 

2.2 

CS31062-050 

1.5 

2.00 X 10-3 

0.54 

2.42 X 10* 

0.32 

2.92 X 10* 


96.6 

20 

4.8 

HD 187861 

1.5 

2.00 X 10-3 

0.59 

4.78 X 10* 

0.31 

6.49 X 10* 


16.6 

8 

2.1 

HD224959 

1.5 

2.00 X 10-3 

0.54 

2.42 X 10* 

0.32 

2.96 X 10* 

1.42 X 103 

64.6 

8 

8.1 

HE0131-3953 

1.5 

2.00 X 10-3 

0.61 

6.71 X 10* 

0.28 

9.88 X 10* 


13.2 

4 

3.3 

HE0143-0441 

1.5 

2.00 X 10-3 

0.59 

4.78 X 10* 

0.31 

6.43 X 10* 


37.2 

7 

5.3 

HE0338-3945 

1.5 

6.66 X 10-* 

0.49 

2.45 X 10* 

0.30 

2.88 X 10* 


183.7 

17 

10.8 

HEl 105+0027 

1.5 

2.00 X 10-3 

0.56 

3.40 X 10* 

0.32 

4.28 X 10* 


23.0 

4 

5.8 

HE1305+0007 

1.5 

6.66 X 10-* 

0.54 

2.42 X 10* 

0.32 

2.97 X 10* 


56.0 

10 

5.6 

HE2148-1247 

1.5 

2.00 X 10-3 

0.51 

2.44 X 10* 

0.31 

2.89 X 10* 


204.7 

11 

18.6 

HE2258-6358 

1.5 

6.66 X 10-* 

0.59 

4.78 X 10* 

0.31 

6.62 X 10* 


293.5 

18 

16.3 

LP625-44 

1.7 

4.00 X 10-3 

0.49 

2.34 X 10* 

0.34 

2.86 X 10* 

4.03 X 103 

173.2 

17 

10.2 

SDSSJ0912+0216 

1.5 

2.00 X 10-3 

0.46 

2.47 X 10* 

0.28 

2.87 X 10* 


425.0 

14 

30.4 

CEMP-j/ur stars 

BD+04°2466 

1.1 

2.00 X 10-3 

0.76 

3.19 X 103 

0.09 

4.27 X 103 

4.29 X 103 

13.4 

7 

1.9 

HD 13826 

1.7 

0 

0.79 

3.90 X 103 

0.07 

8.85 X 103 


3.7 

3 

1.2 

HD187216 

1.5 

5.00 X 10-* 

0.56 

1.71 X 10* 

0.32 

2.33 X 10* 


28.7 

4 

7.2 

HD201626 

1.4 

4.00 X 10-3 

0.71 

1.34 X 103 

0.18 

2.33 X 103 

4.75 X 103 

7.5 

5 

1.5 

HD5223 

1.1 

1.50 X 10-3 

0.76 

1.79 X 10* 

0.17 

2.44 X 10* 


15.0 

9 

1.7 

HE0024-2523 

0.9 

2.00 X 10-3 

0.81 

4.19 X 103 

0.12 

4.75 X 103 

3.33 

35.4 

4 

00 

bo 

HE0212-0557 

1.5 

2.00 X 10-3 

0.59 

4.78 X 10* 

0.31 

6.51 X 10* 


67.9 

5 

13.6 

HE0231-4016 

0.9 

4.00 X 10-3 

0.79 

8.42 X 103 

0.07 

8.81 X 103 


1.4 

4 

0.4 

HE0430-4404 

0.9 

4.00 X 10-3 

0.84 

5.88 X 103 

0.09 

6.25 X 103 


1.0 

2 

0.5 

HE1430-1123 

1.5 

2.00 X 10-3 

0.76 

1.83 X 103 

0.08 

3.97 X 103 


4.6 

2 

2.3 

HE 1434-1442 

1.0 

3.00 X 10-3 

0.79 

2.06 X 103 

0.05 

2.73 X 103 


17.7 

3 

5.9 

HE2150-0825 

0.9 

1.50 X 10-3 

0.84 

5.88 X 103 

0.09 

6.25 X 103 


1.1 

5 

0.2 

HE2232-0603 

0.9 

1.00 X 10-3 

0.84 

5.88 X 103 

0.09 

6.25 X 103 


9.7 

4 

2.4 

SDSS0924+40 

1.5 

2.00 X 10-3 

0.51 

2.44 X 10* 

0.31 

2.89 X 10* 


14.0 

2 

7.0 

SDSS2047+00 

1.5 

6.66 X 10-* 

0.54 

1.36 X 10’ 

0.07 

3.05 X 103 


0.7 

3 

0.2 


Notes. Masses are expressed in units of Mg, periods in days. 
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Table A.3. Physical parameters of the model CEMP-s stars with number of degrees of freedom v > 2 computed with model set B 
in which we adopt an efficient BHL model of wind mass transfer and efficient angular momentum loss. 


ID 

Afi,i 

iffpMZ 

M2,i 


AMacc 

Pf 

Pf-Paper I 

xv mm 

y 

y3. /v 

XI min ' 

CEMP-j/«r stars 

CS22880-074 

1.1 

1.00 X 10-" 

0.79 

7.10 X 10* 

0.05 

5.42 X 10* 


16.6 

10 

1.7 

CS22942-019 

1.8 

2.00 X 10-“* 

0.56 

1.59 X 10* 

0.32 

3.48 X 103 

3.23 X 103 

21.0 

9 

2.3 

CS22964-161A 

1.6 

2.00 X 10^3 

0.79 

1.30 X 103 

0.04 

1.56 X 103 

1.30 X 103 

19.6 

8 

2.5 

CS22964-161B 

1.6 

2.00 X 10^3 

0.71 

1.30 X 10’ 

0.02 

1.56 X 103 

1.30 X 103 

11.8 

9 

1.3 

CS30301-015 

2.4 

4.00 X 10^3 

0.79 

2.17 X 103 

0.04 

4.09 X 10* 


11.0 

8 

1.4 

HD 196944 

1.2 

1.00 X 10^3 

0.79 

6.92 X 10* 

0.04 

4.63 X 10* 


13.6 

12 

1.1 

HD 198269 

1.5 

3.00 X 10-3 

0.79 

4.56 X 10* 

0.08 

2.29 X 10* 

1.17 X 103 

9.4 

5 

1.9 

HE0202-2204 

0.9 

4.00 X 10-3 

0.79 

1.68 X 103 

0.07 

1.51 X 103 


5.1 

7 

0.7 

HEl 1354-0139 

0.9 

6.66 X 10-* 

0.81 

1.05 X 103 

0.03 

9.47 X 10* 


3.2 

6 

0.5 

HE2158-0348 

1.5 

3.00 X 10-3 

0.74 

2.31 X 10* 

0.15 

1.17 X 10* 


16.5 

7 

2.4 


CEMP-j/r stars 


CS22881-036 

1.5 

6.66 X 10-* 

0.71 

2.93 X 10’ 

0.14 

6.86 X 10’ 


9.4 

6 

1.6 

CS22898-027 

1.1 

1.00 X 10-3 

0.64 

3.30 X 103 

0.22 

2.24 X 103 


114.0 

11 

10.4 

CS22948-027 

1.5 

2.00 X 10-3 

0.61 

5.98 X 103 

0.25 

2.56 X 103 

3.69 X 103 

29.5 

9 

3.3 

CS29497-030 

1.6 

2.00 X 10-3 

0.51 

1.19 X 10* 

0.30 

2.31 X 103 

3.00 X 103 

145.8 

16 

9.1 

CS29526-110 

1.6 

2.00 X 10-3 

0.54 

1.19 X 10* 

0.30 

3.07 X 103 


25.2 

6 

4.2 

CS31062-012 

1.5 

2.00 X 10-3 

0.54 

8.60 X 103 

0.23 

4.69 X 103 


17.6 

8 

2.2 

CS31062-050 

1.6 

2.00 X 10-3 

0.54 

1.19 X 10* 

0.30 

3.07 X 103 


115.9 

20 

5.8 

HD 187861 

1.5 

2.00 X 10-3 

0.59 

1.20 X 10* 

0.26 

4.35 X 103 


17.3 

8 

2.2 

HD224959 

1.8 

4.00 X 10-3 

0.54 

1.60 X 10* 

0.31 

2.59 X 103 

1.41 X 103 

70.5 

8 

8.8 

HE0131-3953 

1.5 

2.00 X 10-3 

0.61 

5.98 X 103 

0.25 

2.59 X 103 


13.6 

4 

3.4 

HE0143-0441 

1.5 

2.00 X 10-3 

0.54 

1.21 X 10* 

0.24 

3.65 X 103 


39.5 

7 

5.6 

HE0338-3945 

1.1 

1.50 X 10-3 

0.59 

3.35 X 103 

0.22 

2.05 X 103 


211.8 

17 

12.5 

HE 1105+0027 

1.5 

2.00 X 10-3 

0.59 

1.20 X 10* 

0.26 

4.38 X 103 


24.2 

4 

6.0 

HE1305+0007 

1.1 

1.00 X 10-3 

0.66 

3.28 X 103 

0.22 

2.30 X 103 


64.2 

10 

6.4 

HE2148-1247 

1.6 

2.00 X 10-3 

0.54 

1.19 X 10* 

0.30 

3.07 X 103 


243.8 

11 

22.2 

HE2258-6358 

1.5 

1.00 X 10-3 

0.59 

1.20 X 10* 

0.26 

4.28 X 103 


301.9 

18 

16.8 

LP625-44 

1.8 

4.00 X 10-3 

0.54 

1.60 X 10* 

0.31 

2.61 X 103 

3.65 X 103 

188.9 

17 

11.1 

SDSSJ0912+0216 

1.5 

2.00 X 10-3 

0.54 

8.60 X 103 

0.23 

4.69 X 103 


454.2 

14 

32.4 

CEMP-j/ur stars 

BD+04°2466 

0.9 

2.00 X 10-3 

0.76 

8.49 X 103 

0.12 

6.64 X 103 

4.25 X 103 

13.5 

7 

1.9 

HD 13826 

1.7 

0 

0.76 

8.78 X 10* 

0.05 

2.95 X 10* 


4.0 

3 

1.3 

HD187216 

1.1 

2.00 X 10-* 

0.66 

3.28 X 103 

0.22 

2.07 X 103 


23.8 

4 

5.9 

HD201626 

1.4 

4.00 X 10-3 

0.66 

1.71 X 10* 

0.16 

8.35 X 103 

3.55 X 103 

7.8 

5 

1.6 

HD5223 

1.1 

1.50 X 10-3 

0.74 

1.81 X 10* 

0.12 

1.29 X 10* 


14.6 

9 

1.6 

HE0024-2523 

0.9 

4.00 X 10-3 

0.76 

6.01 X 103 

0.16 

5.23 X 103 

3.48 

37.0 

4 

9.3 

HE0212-0557 

1.5 

2.00 X 10-3 

0.59 

1.20 X 10* 

0.26 

4.34 X 103 


70.1 

5 

14.0 

HE0231-4016 

0.9 

1.50 X 10-3 

0.76 

8.49 X 103 

0.12 

6.80 X 103 


1.1 

4 

0.3 

HE0430-4404 

0.9 

3.00 X 10-3 

0.74 

8.55 X 103 

0.11 

6.73 X 103 


0.7 

2 

0.4 

HE1430-1123 

1.0 

4.00 X 10-3 

0.81 

2.57 X 10* 

0.09 

2.31 X 10* 


4.5 

2 

2.2 

HE1434-1442 

1.0 

3.00 X 10-3 

0.84 

7.19 X 10* 

0.05 

6.34 X 10* 


18.3 

3 

6.1 

HE2150-0825 

0.9 

4.00 X 10-3 

0.79 

4.22 X 103 

0.20 

4.08 X 103 


0.8 

5 

0.2 

HE2232-0603 

0.9 

6.66 X 10-* 

0.76 

6.01 X 103 

0.16 

5.23 X 103 


8.8 

4 

2.2 

SDSS0924+40 

1.6 

2.00 X 10-3 

0.51 

1.19 X 10* 

0.30 

2.31 X 103 


16.0 

2 

8.0 

SDSS2047+00 

1.1 

1.50 X 10-3 

0.51 

3.85 X 10* 

0.05 

2.04 X 10* 


0.8 

3 

0.3 


Notes. Masses are expressed in units of Mq, periods in days. 
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Table A.4. Physical parameters of the model CEMP-s stars with number of degrees of freedom v > 2 computed with model set C 
in which we adopt the WRLOF wind-accretion efficiency and efficient angular momentum loss. 


ID 

Ml-, 

iffpMZ 


Mz,, 

J 

P- 


AA/acc 

J 

Pf 


- 

^ min 

V 

. /y 

^ mm • 

CEMP-i/nr stars 

CS22880-074 

1.1 

1.00 

X 

10 


0.79 

6.33 

X 

10^ 

0.09 

4.94 

X 

10® 

16.8 

10 

1.7 

CS22942-019 

1.8 

2.00 

X 

10 

-4 

0.54 

6.38 

X 

lO'* 

0.29 

1.47 

X 

10® 

20.9 

9 

2.3 

CS22964-161A 

1.6 

2.00 

X 

10 

-3 

0.79 

2.60 

X 

10® 

0.04 

3.05 

X 

10® 

19.7 

8 

2.5 

CS22964-161B 

1.6 

2.00 

X 

10 

-3 

0.71 

2.60 

X 

10® 

0.03 

3.05 

X 

10® 

11.5 

9 

1.3 

CS30301-015 

2.4 

2.00 

X 

10 

-3 

0.81 

1.37 

X 

lO'* 

0.05 

2.28 

X 

10® 

11.0 

8 

1.4 

HD 196944 

1.2 

1.00 

X 

10 

-3 

0.79 

8.71 

X 

10^ 

0.07 

5.89 

X 

10® 

13.9 

12 

1.2 

HD198269 

1.0 

2.00 

X 

10 

-3 

0.84 

2.03 

X 

10® 

0.07 

1.77 

X 

10® 

9.2 

5 

1.8 

HE0202-2204 

0.9 

2.00 

X 

10 

-3 

0.81 

1.67 

X 

10® 

0.07 

1.51 

X 

10® 

5.1 

7 

0.7 

HE1135-H0139 

0.9 

1.00 

X 

10 

-3 

0.79 

4.74 

X 

10® 

0.09 

4.43 

X 

10® 

2.9 

6 

0.5 

HE2158-0348 

1.5 

3.00 

X 

10 

-3 

0.71 

3.69 

X 

10® 

0.15 

1.76 

X 

10® 

15.4 

7 

2.2 


CEMP-i/r stars 


CS22881-036 

1.5 

6.66 

X 

10 


0.71 

3.69 

X 

10® 

0.15 

1.77 

X 

10® 

10.1 

6 

1.7 

CS22898-027 

1.5 

6.66 

X 

10 

-4 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

101.6 

11 

9.2 

CS22948-027 

1.5 

2.00 

X 

10 

-3 

0.64 

1.88 

X 

10® 

0.24 

8.29 

X 

10® 

29.6 

9 

3.3 

CS29497-030 

1.5 

2.00 

X 

10 

-3 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

129.6 

16 

8.1 

CS29526-110 

1.5 

2.00 

X 

10 

-3 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

22.5 

6 

3.8 

CS31062-012 

1.5 

2.00 

X 

10 

-3 

0.59 

1.35 

X 

10® 

0.26 

5.39 

X 

10® 

17.5 

8 

2.2 

CS31062-050 

1.5 

2.00 

X 

10 

-3 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

98.3 

20 

4.9 

HD187861 

1.5 

2.00 

X 

10 

-3 

0.56 

9.59 

X 

10® 

0.28 

3.63 

X 

10® 

16.6 

8 

2.1 

HD224959 

1.5 

2.00 

X 

10 

-3 

0.54 

4.84 

X 

10® 

0.31 

1.74 

X 

10® 

65.6 

8 

8.2 

HE0131-3953 

1.5 

2.00 

X 

10 

-3 

0.61 

1.34 

X 

10® 

0.27 

5.79 

X 

10® 

13.3 

4 

3.3 

HE0143-0441 

1.5 

2.00 

X 

10 

-3 

0.56 

9.59 

X 

10® 

0.28 

3.65 

X 

10® 

37.2 

7 

5.3 

HE0338-3945 

1.5 

6.66 

X 

10 

-4 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

188.6 

17 

11.1 

HE 1105+0027 

1.5 

2.00 

X 

10 

-3 

0.56 

6.79 

X 

10® 

0.31 

2.66 

X 

10® 

23.0 

4 

5.8 

HE 1305+0007 

1.5 

6.66 

X 

10 

-4 

0.54 

4.84 

X 

10® 

0.31 

1.73 

X 

10® 

56.7 

10 

5.7 

HE2148-1247 

1.5 

2.00 

X 

10 

-3 

0.56 

4.81 

X 

10® 

0.32 

1.92 

X 

10® 

208.9 

11 

19.0 

HE2258-6358 

1.5 

6.66 

X 

10 

-4 

0.59 

6.75 

X 

10® 

0.32 

2.81 

X 

10® 

293.4 

18 

16.3 

LP625-44 

1.7 

4.00 

X 

10 

-3 

0.54 

6.52 

X 

10® 

0.35 

1.86 

X 

10® 

177.6 

17 

10.4 

SDSS10912+0216 

1.5 

2.00 

X 

10 

-3 

0.49 

6.92 

X 

10® 

0.27 

2.03 

X 

10® 

429.1 

14 

30.7 

CEMP-i/ur stars 

BD+04°2466 

1.1 

1.50 

X 

10 


0.74 

4.54 

X 

10® 

0.10 

3.23 

X 

10® 

13.2 

7 

1.9 

HD13826 

1.7 


0 



0.76 

1.11 

X 

10® 

0.06 

3.38 

X 

10® 

3.8 

3 

1.3 

HD187216 

1.5 

5.00 

X 

10 

-4 

0.56 

4.81 

X 

10® 

0.32 

1.70 

X 

10® 

28.7 

4 

7.2 

HD201626 

1.4 

4.00 

X 

10 

-3 

0.71 

2.67 

X 

10® 

0.18 

1.49 

X 

10® 

7.5 

5 

1.5 

HD5223 

1.1 

1.50 

X 

10 

-3 

0.76 

2.54 

X 

10® 

0.16 

1.99 

X 

10® 

15.1 

9 

1.7 

HE0024-2523 

1.1 


0 



0.69 

7.29 

X 

10® 

0.14 

5.39 

X 

10® 

47.7 

4 

11.9 

HE0212-0557 

1.5 

2.00 

X 

10 

-3 

0.59 

9.53 

X 

10® 

0.29 

3.91 

X 

10® 

67.9 

5 

13.6 

HE0231-4016 

0.9 

1.50 

X 

10 

-3 

0.79 

8.42 

X 

10® 

0.08 

6.62 

X 

10® 

1.1 

4 

0.3 

HE0430-4404 

0.9 

4.00 

X 

10 

-3 

0.79 

8.42 

X 

10® 

0.08 

6.62 

X 

10® 

0.7 

2 

0.3 

HE1430-1123 

1.5 

2.00 

X 

10 

-3 

0.74 

5.18 

X 

10® 

0.09 

2.45 

X 

10® 

4.5 

2 

2.2 

HE1434-1442 

1.0 

4.00 

X 

10 

-3 

0.79 

2.90 

X 

10® 

0.05 

2.46 

X 

10® 

17.8 

3 

5.9 

HE2150-0825 

0.9 

3.00 

X 

10 

-3 

0.84 

4.16 

X 

10® 

0.14 

3.95 

X 

10® 

2.2 

5 

0.4 

HE2232-0603 

0.9 

6.66 

X 

10 

-4 

0.84 

5.88 

X 

10® 

0.11 

5.20 

X 

10® 

10.1 

4 

2.5 

SDSS0924+40 

1.5 

2.00 

X 

10 

-3 

0.54 

4.84 

X 

10® 

0.31 

1.75 

X 

10® 

14.4 

2 

7.2 

SDSS2047+00 

1.5 

6.66 

X 

10 

-4 

0.61 

5.33 

X 

10® 

0.08 

1.90 

X 

10® 

0.7 

3 

0.2 


Notes. Masses are expressed in units of M©, periods in days. 
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Table A.5. Physical parameters of the model CEMP-s stars with number of degrees of freedom v < 1 computed with model set A. 


ID 

Ml-, 

Mvml 

M2,i 

Pi 

AA/acc 

Pf 

Pf-Paper I 

- 

^ min 

V 

1 , 
y~ • /v 

mm ' 

CEMP-s/r stars 

BS16080-175 

2.9 

4.00 X 10-' 

0.59 

9.29 X 103 

0.26 

2.55 X 10'' 


20.7 

1 

20.7 

BS17436-058 

1.1 

4.00 X 10-3 

0.74 

3.61 X lO'* 

0.15 

4.73 X lO"* 


18.4 

1 

18.4 

CEMP-s/ur stars 

CS22945-017 

2.0 

2.00 X IQ-'* 

0.74 

1.66 X 103 

0.13 

4.50 X 103 


3.7 

-1 


CS22956-028 

1.3 

0 

0.64 

1.25 X lO'* 

0.20 

1.77 X lO"* 

1.37 X 103 

14.8 

1 

14.8 

CS29509-027 

2.9 

1.50 X 10-3 

0.69 

3.25 X 103 

0.11 

5.28 X 103 

1.65 X 103 

0.8 

-1 


CS30338-089 

1.5 

6.66 X lO-'* 

0.59 

4.78 X lO'* 

0.31 

6.53 X lO-* 


6.1 

-1 


HE0012-1441 

2.6 

4.00 X 10-3 

0.76 

4.74 X 103 

0.06 

1.58 X lO"* 


0.2 

-1 


HE0206-1916 

1.5 

6.66 X lO-'* 

0.64 

9.42 X lO'* 

0.24 

1.57 X 103 


2.5 

-1 


HE0322-1504 

1.8 

4.00 X 10-3 

0.49 

2.29 X lO'* 

0.36 

3.13 X lO"* 


1.7 

-1 


HE0441-0652 

1.2 

3.00 X 10-3 

0.79 

1.55 X 103 

0.05 

1.90 X 103 


0.2 

-1 


HE0507-1430 

6.0 

0 

0.56 

6.77 X 103 

0.26 

3.83 X lO-* 

5.16 X 103 

0.9 

-1 


HE0534-4548 

2.3 

1.00 X 10-3 

0.81 

2.20 X 103 

0.02 

9.10 X 103 


0.0 

-2 


HE1045-1434 

6.0 

0 

0.49 

1.36 X lO'* 

0.37 

5.56 X lO-* 


17.1 

-1 


HEl 157-0518 

1.5 

6.66 X lO-'* 

0.59 

4.78 X lO'* 

0.31 

6.55 X lO-* 


1.7 

-1 


HE1429-0551 

1.3 

1.00 X 10-3 

0.71 

1.37 X 103 

0.17 

2.29 X 103 


0.5 

-1 


HE1443-H0113 

1.7 

0 

0.74 

2.79 X 103 

0.09 

4.67 X 103 


0.1 

-1 


HE1447-H0102 

1.7 

6.66 X lO-'* 

0.49 

2.34 X lO'* 

0.34 

2.94 X lO"* 


73.7 

-1 


HE1523-1155 

1.7 

4.00 X 10-3 

0.76 

1.75 X 103 

0.11 

4.35 X 103 


1.2 

-1 


HE1528-0409 

1.5 

1.50 X 10-3 

0.54 

2.42 X lO'* 

0.32 

2.99 X lO"* 


12.2 

-1 


HE2221-0453 

1.8 

4.00 X 10-3 

0.74 

1.73 X 103 

0.11 

5.33 X 103 


1.3 

-1 


HE2227-4044 

1.0 

3.00 X 10-3 

0.81 

4.57 X 103 

0.10 

5.55 X 103 


0.3 

1 

0.3 

HE2228-0706 

1.7 

2.00 X 10-3 

0.51 

2.33 X lO'* 

0.36 

2.87 X lO"* 


7.6 

-1 


HE2240-0412 

1.0 

3.00 X 10-3 

0.79 

6.50 X 103 

0.09 

7.96 X 103 


1.2 

-1 


SDSS1707+58 

1.8 

4.00 X 10-3 

0.49 

2.29 X lO'* 

0.36 

2.87 X lO"* 


22.3 

1 

22.3 


Notes. Masses are expressed in units of M©, periods in days. 

Table A.6. Physical parameters of the model CEMP-s stars with number of degrees of freedom v < 1 computed with model set B. 


ID 

Mi,i 

MpMZ 

M2,i 

Pi 


Pi 

Pf-Paper I 

- 

^ min 

y 

y3 . /y 
/V min' 

CEMP-s/r stars 

BS 16080-175 

3.0 

0 

0.59 

3.64 X 10" 

0.23 

2.86 X 103 


20.8 

1 

20.8 

BS 17436-058 

1.5 

1.50 X 10-3 

0.74 

2.91 X 103 

0.15 

7.30 X 10* 


18.4 

1 

18.4 

CEMP-s/ur stars 

CS22945-017 

2.0 

2.00 X 10-" 

0.71 

3.74 X 10' 

0.16 

3.65 X 10' 


3.4 

-1 


CS22956-028 

1.3 

0 

0.69 

1.23 X 10" 

0.19 

7.36 X 103 

1.05 X 103 

14.3 

1 

14.3 

CS29509-027 

3.0 

0 

0.66 

6.41 X 103 

0.15 

2.08 X 10' 

1.67 X 103 

0.4 

-1 


CS30338-089 

1.1 

1.50 X 10-3 

0.66 

4.63 X 103 

0.22 

3.35 X 103 


5.8 

-1 


HE0012-1441 

2.6 

4.00 X 10-3 

0.76 

2.12 X 103 

0.05 

3.38 X 10" 


0.3 

-1 


HE0206-1916 

1.1 

1.50 X 10-3 

0.74 

6.41 X 103 

0.21 

5.29 X 103 


1.9 

-1 


HE0322-1504 

1.8 

4.00 X 10-3 

0.54 

1.60 X 10" 

0.31 

2.40 X 103 


2.3 

-1 


HE0441-0652 

1.5 

3.00 X 10-3 

0.79 

1.29 X 103 

0.03 

6.07 X 10" 


0.1 

-1 


HE0507-1430 

6.0 

0 

0.59 

2.14 X 103 

0.27 

2.51 X 103 

3.58 X 103 

1.0 

-1 


HE0534-4548 

2.8 

6.66 X 10-" 

0.81 

5.76 X 103 

0.02 

8.47 X 10" 


0.0 

-2 


HE1045-1434 

6.0 

0 

0.56 

3.02 X 103 

0.29 

2.61 X 103 


20.9 

-1 


HE1157-0518 

1.1 

1.50 X 10-3 

0.66 

3.28 X 103 

0.22 

2.30 X 103 


1.7 

-1 


HE1429-0551 

1.3 

1.00 X 10-3 

0.66 

1.24 X 10" 

0.18 

6.78 X 103 


0.5 

-1 


HE1443+0113 

1.3 

6.66 X 10-" 

0.74 

2.42 X 10" 

0.11 

1.45 X 10" 


0.3 

-1 


HE1447+0102 

1.2 

2.00 X 10-3 

0.61 

3.24 X 103 

0.35 

1.39 X 103 


79.1 

-1 


HE1523-1155 

1.7 

4.00 X 10-3 

0.76 

2.78 X 103 

0.12 

3.76 X 10' 


1.2 

-1 


HE1528-0409 

1.6 

1.00 X 10-3 

0.54 

1.19 X 10" 

0.30 

3.01 X 103 


13.1 

-1 


HE2221-0453 

1.8 

4.00 X 10-3 

0.74 

4.33 X 10" 

0.12 

1.34 X 10" 


1.3 

-1 


HE2227-4044 

1.0 

1.50 X 10-3 

0.81 

7.24 X 10" 

0.05 

6.26 X 10" 


0.2 

1 

0.2 

HE2228-0706 

1.7 

3.00 X 10-3 

0.59 

1.62 X 10" 

0.26 

4.49 X 103 


9.7 

-1 


HE2240-0412 

1.5 

2.00 X 10-3 

0.74 

2.06 X 103 

0.04 

3.67 X 10' 


0.8 

-1 


SDSS1707+58 

1.8 

4.00 X 10-3 

0.54 

1.60 X 10" 

0.31 

2.61 X 103 


24.9 

1 

24.9 


Notes. Masses are expressed in units of M©, periods in days. 
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Table A.7. Physical parameters of model CEMP-i with number of degrees of freedom v < 1 computed with model set C 


ID 

iWi,i 

Mfuz 


M2, 


P, 


AMacc 

i 



min 

V 

- /v 

^ min' 

CEMP-i/r stars 

BS16080-175 

2.9 

1.00 

X 

10 


0.56 

2.09 

X 

10^ 

0.26 

1.57 

X 

10^ 

20.7 

1 

20.7 

BS17436-058 

1.1 

4.00 

X 

10- 

-3 

0.71 

2.57 

X 

lO-* 

0.16 

2.00 

X 

10^ 

18.4 

1 

18.4 

CEMP-i/ur stars 

CS22945-017 

2.0 

2.00 

X 

10- 


0.69 

2.98 

X 

lO'^ 

0.15 

7.49 

X 

10’ 

4.0 

-1 


CS22956-028 

1.3 


0 



0.64 

2.49 

X 

10-* 

0.19 

1.37 

X 

10^ 

14.8 

1 

14.8 

CS29509-027 

3.0 


0 



0.69 

3.60 

X 

lo-* 

0.11 

5.09 

X 

10^ 

0.5 

-1 


CS30338-089 

1.5 

6.66 

X 

10 

-4 

0.59 

9.53 

X 

lo-* 

0.29 

3.90 

X 

10^ 

6.1 

-1 


HE0012-1441 

2.6 

4.00 

X 

10 

-3 

0.74 

1.89 

X 

10-* 

0.07 

2.55 

X 

102 

0.3 

-1 


HE0206-1916 

1.5 

1.00 

X 

10 

-3 

0.66 

2.35 

X 

lo-* 

0.18 

1.04 

X 

10^ 

2.4 

-1 


HE0322-1504 

1.8 

4.00 

X 

10- 

-3 

0.54 

9.01 

X 

lo-* 

0.36 

2.09 

X 

10^ 

1.9 

-1 


HE0441-0652 

1.5 

3.00 

X 

10 

-3 

0.81 

2.86 

X 

10^ 

0.04 

7.12 

X 

10* 

0.1 

-1 


HE0507-1430 

4.0 


0 



0.54 

1.02 

X 

10® 

0.32 

2.76 

X 

10^ 

1.1 

-1 


HE0534-4548 

1.9 

6.66 

X 

10 

-4 

0.84 

1.32 

X 

10® 

0.02 

4.51 

X 

10^ 

0.0 

-2 


HE1045-1434 

2.0 

4.00 

X 

10- 

-3 

0.51 

1.23 

X 

10^ 

0.35 

2.10 

X 

10^ 

24.6 

-1 


HE 1157-0518 

1.5 

6.66 

X 

10 

-4 

0.59 

6.75 

X 

lo-* 

0.32 

2.83 

X 

10^ 

1.7 

-1 


HE 1429-0551 

1.3 

1.00 

X 

10 

-3 

0.69 

2.75 

X 

10^ 

0.15 

1.55 

X 

10^ 

0.7 

-1 


HE1443-r0113 

1.7 


0 



0.74 

7.01 

X 

10^ 

0.09 

2.60 

X 

10^ 

0.1 

-1 


HE1447+0102 

1.7 

6.66 

X 

10 

-4 

0.54 

6.52 

X 

lo-* 

0.35 

1.82 

X 

10^ 

74.5 

-1 


HE1523-1155 

1.5 

4.00 

X 

10 

-3 

0.74 

3.67 

X 

10^ 

0.15 

1.81 

X 

10^ 

1.3 

-1 


HE1528-0409 

1.5 

1.50 

X 

10 

-3 

0.54 

4.84 

X 

10-* 

0.31 

1.72 

X 

10^ 

12.3 

-1 


HE2221-0453 

1.8 

4.00 

X 

10 

-3 

0.76 

6.84 

X 

10^ 

0.12 

2.22 

X 

10^ 

1.3 

-1 


HE2227-4044 

1.5 

1.00 

X 

10 

-3 

0.76 

5.78 

X 

10^ 

0.06 

3.89 

X 

102 

0.2 

1 

0.2 

HE2228-0706 

1.7 

2.00 

X 

10 

-3 

0.54 

6.52 

X 

lo-* 

0.35 

1.85 

X 

10^ 

7.8 

-1 


HE2240-0412 

1.5 

2.00 

X 

10 

-3 

0.76 

4.09 

X 

10^ 

0.05 

1.51 

X 

102 

1.0 

-1 


SDSS1707-H58 

1.8 

4.00 

X 

10 

-3 

0.51 

9.06 

X 

10-* 

0.34 

2.04 

X 

10^ 

23.2 

1 

23.2 


Notes. Masses are expressed in units of Mq, periods in days. 
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Table A.8. Confidence ranges of the input parameters of the modelled CEMP-i stars with;^f^/v < 3 and v >2 computed with model 
set A with WRLOF wind-accretion efficiency and spherically symmetric wind. 


ID 




iffpMz/10 

-3 





P, 



best 

min 

max 

best 

min 

max 

best 

min 

max 

best 

min 

max 

CEMP-j/nr 





CS22880-074 

1.10 

0.85 

1.15 

1.00 

0.00 

4.00 

0.74 

0.73 

0.85 

4.54 X 10" 

6.64 X 10" 

1.64 X lO'* 



1.45 

1.55 








1.56 X 10" 

3.47 X 10" 

CS22942-019 

1.80 

1.75 

1.85 

0.20 

0.00 

0.58 

0.54 

0.48 

0.68 

6.48 X 10“* 

8.73 X 10" 

1.70 X 10" 

CS22964-161A 

1.60 

0.95 

1.65 

2.00 

0.58 

4.00 

0.79 

0.75 

0.85 

2.74 X 10^ 

2.18 X 10" 

4.87 X 10" 

CS22964-161B 

1.60 

0.95 

1.65 

2.00 

0.58 

4.00 

0.71 

0.68 

0.85 

2.74 X 10^ 

2.18 X 10" 

4.87 X 10" 

CS30301-015 

2.50 

1.95 

2.75 

4.00 

0.58 

4.00 

0.76 

0.70 

0.85 

5.14 X 10" 

2.90 X 10" 

6.47 X 10" 



1.15 

1.25 








1.30 X 10" 

2.27 X 10" 

HD 196944 

1.40 

0.95 

1.65 

1.00 

0.30 

4.00 

0.81 

0.75 

0.88 

5.84 X 10" 

1.38 X 10" 

9.92 X 10" 

HD198269 

1.00 

0.95 

1.65 

4.00 

0.80 

4.00 

0.79 

0.73 

0.85 

1.64 X 10" 

1.31 X 10" 

2.05 X 10" 



1.95 

2.15 










HE0202-2204 

0.90 

0.85 

1.15 

3.00 

0.00 

4.00 

0.79 

0.75 

0.88 

1.15 X 10" 

4.67 X 10" 

1.05 X 10“* 



1.45 

1.55 








2.06 X 10" 

2.60 X 10" 

HE1135-r0139 

0.90 

0.85 

1.95 

1.00 

0.00 

4.00 

0.79 

0.75 

0.88 

3.35 X 10" 

4.67 X 10" 

7.29 X 10" 

HE2158-0348 

1.50 

1.15 

1.65 

3.00 

1.75 

4.00 

0.69 

0.65 

0.75 

1.18 X 10" 

5.79 X 10" 

9.17 X 10" 












9.34 X 10^ 

1.48 X 10" 

CEMP-i/r 





CS22881-036 

1.50 

1.45 

1.55 

0.67 

0.05 

2.50 

0.66 

0.65 

0.80 

1.35 X 10" 

6.25 X 10" 

9.90 X 10" 



1.05 

1.15 








1.38 X 10^ 

1.71 X 10" 

CS31062-012 

1.50 

1.35 

1.65 

2.00 

0.58 

3.50 

0.56 

0.45 

0.65 

6.92 X 10“* 

9.65 X 10" 

1.21 X 10" 

HD187861 

1.50 

1.35 

1.65 

2.00 

0.58 

3.50 

0.59 

0.50 

0.65 

4.90 X 10“* 

9.59 X 10" 

1.21 X 10" 

CEMP-i/ur 





BD+04°2466 

1.10 

0.85 

1.15 

2.00 

0.00 

4.00 

0.76 

0.73 

0.85 

3.19 X 10" 

5.70 X 10" 

5.27 X 10" 



1.45 

1.65 








9.90 X 10" 

2.37 X 10" 

HD13826 

1.70 

1.45 

1.85 

0.00 

0.00 

0.30 

0.79 

0.68 

0.88 

3.90 X 10" 

4.70 X 10" 

1.16 X 10“* 



0.85 

0.95 








1.62 X 10" 

2.57 X 10" 

HD201626 

1.40 

1.05 

2.15 

4.00 

0.80 

4.00 

0.71 

0.60 

0.78 

1.18 X 10" 

2.31 X 10" 

1.63 X 10" 

HD5223 

1.10 

1.05 

1.15 

1.50 

0.30 

3.50 

0.76 

0.68 

0.83 

1.79 X 10“* 

1.31 X 10" 

2.07 X 10" 



1.45 

1.55 








6.63 X 10" 

1.05 X 10“* 












1.38 X 10^ 

1.88 X 10" 

HE0231-4016 

0.90 

0.85 

1.25 

4.00 

0.05 

4.00 

0.79 

0.73 

0.85 

8.42 X 10" 

4.70 X 10" 

7.44 X 10" 



1.45 

1.65 








9.65 X 10" 

1.53 X 10“* 












1.49 X 10" 

2.27 X 10" 

HE0430-4404 

0.90 

0.85 

1.15 

4.00 

0.00 

4.00 

0.84 

0.70 

0.88 

5.88 X 10" 

4.70 X 10" 

5.30 X 10" 



1.45 

1.55 








1.04 X 10^ 

1.66 X 10-* 












1.49 X 10" 

2.37 X 10" 

HE1430-1123 

1.50 

1.35 

1.65 

2.00 

0.05 

4.00 

0.76 

0.68 

0.85 

1.83 X 10" 

7.26 X 10" 

2.05 X 10" 



0.85 

1.25 








1.24 X 10" 

3.17 X 10" 

HE2150-0825 

0.90 

0.85 

1.15 

1.50 

0.05 

4.00 

0.84 

0.73 

0.85 

5.88 X 10" 

3.33 X 10" 

3.75 X 10" 



1.45 

1.55 








9.65 X 10" 

2.37 X 10" 

HE2232-0603 

0.90 

0.85 

1.15 

1.00 

0.00 

1.75 

0.84 

0.73 

0.85 

5.88 X 10" 

3.33 X 10" 

3.24 X 10" 

SDSS2047-(-00 

1.50 

0.85 

1.75 

0.67 

0.05 

4.00 

0.54 

0.43 

0.83 

1.36 X 10" 

3.33 X 10" 

5.28 X 10" 












9.82 X 10" 

2.41 X 10" 


Notes. Masses are expressed in units of Mq, periods in days. 



